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ABSTRACT 

It has been suggested at various times that the angular diameter of a star may be 
deduced from a study of its light-curve while it is undergoing occultation by the moon; 
but Eddington has shown that observation of the geometrical occultation is almost 
hopelessly complicated by diffraction effects, except for stars of very large angular 
diameter. The present paper points out that the diameter may be deduced from the 
diffraction effect itself, and the question is discussed quantitatively. 

MacMahon has proposed? that the diameters of stars be investi- 
gated by determining the time interval required for their occultation 
by the moon. Eddington drew attention to the diffraction effect 
and estimated that this would cause a point star to have a spurious 
diameter of 07008; consequently, the approach to the problem from 
the viewpoint of geometrical optics is rather hopeless except in the 
case of a star with a very large angular diameter. The observation 
in the latter case has been attempted‘ recently at Cambridge and at 
Paris, the star being Antares. 

It has occurred to the writer that the diffraction effect, which 
seems destined to close the direct line of investigation, may be used 
to forward the inquiry; at least this possibility emerges when the 

' Presented at the sixtieth meeting of the American Astronomical Society, at Ann 
Arbor, Michigan. 

2M.N., 69, 126, 1908-1909. 

3M.N., 69, 178, 1908-1909. 4 Nautical Almanac, 1939, Suppl., p. 1, n. 
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problem is idealized in the following manner: if we assume that the 
effect of the moon, so far as diffraction is concerned, is the same as 
that of a plane screen perpendicular to the line of sight and situated 
so that its edge is tangent to the moon’s limb at the point where oc- 
cultation occurs, the phenomenon is reduced to the familiar problem 
of diffraction at a straight edge. For a point source of light, the in- 
tensity of illumination increases from zero within the geometrical 
shadow of the screen, attains one-fourth the intensity of the free 
wave at the edge, and beyond the shadow exhibits a series of oscil- 
lations which converge rapidly to the intensity of the free wave. 
The point to be raised is simply that the character of the oscillations 
for light from an extended source is not the same as for a point source 
and that it depends, in a predictable manner, on the diameter of the 
source. Actually the problem may be inverted so that if the varia- 
tion of intensity in the fringes is observed, we may estimate the 
diameter of the star. 
For a point source the intensity J is given by® 


SJ = (1 + 2F,)* + (1 + 2F.)*, (1) 
the upper signs corresponding to points outside the shadow, the 
lower to points within it; the unit of intensity is that of the free wave. 
The functions F, and F, are Fresnel’s integrals, 


cos — vdv and sin 
° °o 


which have been computed and tabulated. For the case under dis- 
cussion the scale of the phenomenon is remarkable, as may be seen 
from the following list of the abscissas x (in meters), measured from 
the geometric edge of the shadow, corresponding to the first few 
maxima and minima: 


nila 
N13 


Maxima....... 11.0 $%.4 27.8 
PETE ooo ss. 16.9 24.7 


5 Cf. Preston, Theory of Light (4th ed.), p. 296. 
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The values for the moon’s distance D and the wave length of light 
used in obtaining them were 3.8 X 10° and 4.3 X 10 7 meters, the 


formula® being 


The values of v at the turning-points are given by Fresnel. 

The positions of the fringes predicted by theory have been verified 
by observation many times, but the relative intensities have not 
been substantiated quantitatively until Lyman’s recent work.’ 

The writer has sought to obtain by approximate methods the 
intensity-curves corresponding to circular light sources of finite but 
angularly small diameters. A table of the values of J in equation (1) 
was prepared for values of the argument v differing by Av = o.1. 

The desired effect at the point v = 7) was then found from 


I(%, ”) = Cf (af V, VU, n)dv, (2) 


where 





the limits of integration being 1, + nAv and v7. — nAv; C determines 
the unity of intensity. The integration was performed numerically 
by Simpson’s method for several small values of n. Av = 0.1 cor- 
responds to o’0005 displacement of a point source, so integration 
over the range 2nAv represents a circular source of diameter / 1000 
seconds of arc. In Figure 1 the curves are shown for a point source 
and for one of diameter 07004, the range in v being from o to 3. 
Curves for other values of m have been omitted so as to simplify 
the diagram. 

It will be observed that the two light-curves differ considerably, 
which does not seem, at first sight, to be consistent with Lyman’s® 


° Fresnel, GEuvres completes, p. 323, 1866. 


7 Proc. Nat. Acad. Sci., 16, 71, 1930. 3 Op. cit. 
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verification of classical theory. But while Lyman’s light source was 
1°37 when he confirmed the point-source theory, the other constants 
of his apparatus were such as to cause integration over a range in 
v of 0.0175 only. 

If observations of this effect can be secured, the method promises 
to rival the stellar interferometer in power. This last sentence 
brings us squarely to the crux of the matter, namely, the technical 
difficulties. The observations would consist of measuring and re- 
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FIG. 1 


cording the light intensity as the moment of occultation approaches. 
Since the telescope will pass along the fringes at a rate of 1 meter in 
about ofo015, the instruments must be capable of high-speed per- 
formance as well as high sensitivity. However, these criteria are well 
within the limits of time lag and change of reaction to change of 
stimulus associated with photoelectricity, so that the observation is 
conceivable. Perhaps conjunction of the photoelectric cell and the 
cathode-ray oscillograph would make it a reality. 

The present discussion of the subject is intentionally brief and 
crude, since it does not seem worth while to examine it elaborately 
unless observations of this type are actually attempted. In any 
event, theory would be much restricted by the complexity of con- 
ditions on the limb of the moon, although it may be indicated that 
we need be concerned only with a small surface on the moon, since a 
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cone of light with vertex angle o’0006 is only 1 meter across at that 
distance. But there remain such disturbing possibilities as a reflec- 
tion effect. 

Finally, the observation probably should not be attempted with a 
telescope of aperture much less than 1 meter, so that atmospheric 
effects will be integrated out, a suggestion due to R. O. Redman? and 
based upon consideration of the separation of shadow bands. 


I am indebted to Professor H. N. Russell, who mentioned that 
the geometrical problem was beset with diffraction difficulties. 
PRINCETON UNIVERSITY OBSERVATORY 


August 15, 1938 


9 Personal communication. 








PHOTOELECTRIC OBSERVATION OF DIFFRACTION 
AT THE MOON’S LIMB* 


A. E. WHITFORD 


ABSTRACT 

The occultations of 8 Capricorni and vy Aquarii on September 6 and 7 (respectively), 
1938, were observed with the 1oo-inch telescope and a fast photoelectric system. The 
record was obtained by photographing the screen of a cathode-ray tube with a moving 
film camera. An analysis of the unavoidable noise present in amplifiers predicted fluc 
tuations of approximately the magnitude observed, about 4 per cent of the brightness of 
8 Capricorni. 

A calculation of the diffraction effects to be expected was carried out, with allowance 
for the wide range of wave lengths used. The observations were in satisfactory agree- 
ment with the calculation. The stars were not large enough to change the diffraction 
pattern from that due to a point source by an observable amount. Therefore, no in 
formation on stellar diameters was obtained. With stars of angular diameters greater 
than o’005, the method can be used to measure stellar diameters. 

The possibility of determining the diameter of a star by measuring 
the time for it to pass behind the dark limb of the moon has occurred 
independently to several astronomers in the last thirty years, of 
whom MacMahon! seems to have been the first. The problem im- 
mediately resolves itself into two parts: first, the technical question 
of whether it is feasible to obtain observations of the variation in in- 
tensity of the light of a star occurring in a time interval of about 0.01 
sec; and second, the consideration of the influence of diffraction on 
the intensity changes to be expected. The technical question will be 
taken up first. 

The most favorable occultations during 1938 were those of 8 
Capricorni (3"25, composite spectrum Go-+ Ao) and vy Aquarii 
(4'"52 Ko) on September 6 and 7, respectively. The plan of observa- 
tion was to let the light fall on a photoelectric cell at the focus of the 
roo-inch reflector, amplify the variations in a four-stage alternating- 
current amplifier with a voltage gain of about 10°, and apply the 
output to the vertical deflecting plates of a 3-inch cathode-ray oscil- 
loscope. The oscillations of the luminous spot were photographed 

* Contributions from the Mount Wilson Observatory, Carnegie Institution of Washing 


ton, No. 610. 


t M.N., 69, 126, 1900. 
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with a moving-film camera,’ using an f:1.5 lens of 25 mm focus. The 
cathode-ray spot moved only in a vertical direction; the time co- 
ordinate was supplied by the motion of the film. The film speed was 
450 mm/sec. By using a spiral motion of the drum similar to that 
on seismographs, about 75 feet of trace could be recorded on a 3-foot 
piece of 35-mm film. There was a tenfold reduction between 
cathode-ray screen and film. 

This general scheme has been used to record phenomena much 
faster than those here encountered, but always with a considerably 
greater amount of light than that _ lau 
available from a star even witha | 
large telescope. As the intensity ia Pe R ic low 
of the light is reduced, more am- O 





plification is necessary, and final- 


AN 


| 
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ly the signal is no larger than the 

unavoidable background noise L 

present in all amplifiers. nn 
ee 











Figure 1 shows schematically = 
the input circuit of the amplifier. é i - 
: KOS Gp. 8 Fic. 1.—Input circuit of amplifier 
The condenser C, indicated by 
the dotted lines, represents the sum of the capacities of the photocell 
and connections and the dynamic input capacity of the first tube. 
If K is the sensitivity constant of the photocell in amperes per lumen, 
the signal voltage appearing across the input resistor R when a 
luminous flux LZ is stopped or started is 


E, = KLR. (1) 


It is assumed that the signal can be represented with the required 
accuracy, as the sum of sinusoidal waves of frequencies less than 


f. = 1/2eRC. 


On the assumption that noise produced by external disturbances, 
such as stray electric or magnetic fields, vibration, poor connections, 
and the like, has been eliminated, there will remain inherent and 
unavoidable fluctuations arising principally from thermal agitation 
of electricity in the resistor R. In this case, shot effects of grid and 

2 Through the courtesy of Dr. J. A. Anderson this camera was constructed in the 
machine shop of the 200-inch telescope of the California Institute of Technology. 
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photo currents in the input circuit and “tube noise” can be made 
negligible in comparison with thermal noise by proper choice of the 
first tube and its operating potentials. As shown by Nyquist,’ the 
thermal noise voltage is given by 


Nd 


E = ar f R(f)df , (2) 
Oo 


where f is the frequency, R(f) the resistive component of the in- 
put impedance as a function of frequency, 7 the absolute tem- 
perature, and k& the Boltzmann constant. For our case, R(f) = 
R/(1 + 47°R’C?f?). In the amplifier here used, the interstage coup- 
ling condensers were rather large, so that the low-frequency cutoff 
was at about 1 cycle per second. The high-frequency response was 
kept flat to a frequency several times higher than the critical fre- 
quency f, = 1/27RC determined by the capacity shunt in the input 
circuit. Therefore, the amplifier does not appreciably limit the 
spectrum of noise frequencies, and the noise in the input circuit will 
be faithfully reproduced in the output circuit, where the observation 
is made. Hence, it is permissible merely to insert the value of R(/) 
in equation (2) and carry out the integration. The result is‘ 


e=—. (3) 


A numerical calculation can now be made. The total capacity in 
the input circuit was estimated to be 1o micromicrofarads, or 
10% farads. At 300° K the noise voltage then comes out to be 
2X 10 ° RMS volts. The signal voltage for a given light-intensity 
increases directly with R, but the critical frequency f, at which the 
response drops off seriously, varies inversely as RK. Whether a satis- 
factory compromise is possible depends on the amount of light 
available. The .resistance actually used had a value of 2 X 10! 

3 Phys. Rev., 32, 110, 1928. See also J. B. Johnson, ibid., 97; G. L. Pearson, Physics, 
5, 233, 1934; E. A. and A. G. Johnson, Phys. Rev., 50, 171, 1936. 

4 This result was first derived by Hafstad (Phys. Rev., 44, 201, 1935) for a direct- 
current amplifier using an electrometer tube; but, so far as the writer is aware, it has 
not been pointed out that it applies to wide-band alternating-current amplifiers where 
the high-frequency response is limited only by the input circuit. 
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ohms, making the critical frequency f, about 800 cycles. A candle at 
a distance of 1 km is equivalent’ to a star of visual magnitude +0.8. 
Then it may be calculated that 6 Capricorni, visual magnitude 3.25, 
would, after allowing for losses in the atmosphere and three reflec- 
tions from aluminum mirrors, furnish 3 X 10 7 lumens at the focus 
of the 1o00-inch telescope. For radiation from a black body at 
6000° K, the sensitivity constant of a gas-filled cesium oxide photo- 
cell is about 75 X 10 ° amperes per lumen. Then the signal voltage 
E, = 4.5 X 10 4 volts. The signal-to-noise ratio would then be 


Fs _ 4.5 X 104 
(E2)' 2 2 > 4 IO 5 


= = 99 


ie 


which is a reasonably satisfactory value. For v Aquarii, which gives 
a third as much light, the ratio would be only fair. 

The RCA 921 phototube was selected because of its non-micro- 
phonic construction, its low capacity, and its superior high-frequency 
response, compared with other gas-filled cells. The first two tubes of 
the amplifier were RCA 1603's, operated with the screen at 22.5 volts 
and the plate at about 45 volts, the plate supply being from the same 
go-volt battery as that used for the phototube. The 1603’s were 
chosen primarily because of their low microphonic sensitivity, but 
they had a satisfactorily small grid current as well. It was necessary 
to include a triode stage in the part of the amplifier at the end of the 
telescope in order to get the signal into a relatively low-impedance 
circuit before sending it through the 40-foot cable leading from the 
Cassegrain focus to the oscilloscope and its self-contained vertical 
amplifier. Otherwise the capacity of the cable would have shunted 
out the higher frequencies. A 6L5-G tube was used. 

The observations were made successfully on the two nights as 
planned. One observer at the Cassegrain focus of the telescope 
checked the following of the telescope for some minutes and made 
sure that the star image was central in the focal-plane diaphragm up 
to about 40 seconds before immersion, when the guiding prism was 
pulled out and the light admitted to the photocell. The diaphragm 
had a radius of about 1.2 mm, or 6”; the angular distance to the 


’ Russell, Dugan, and Stewart, Astronomy, p. 612, Ginn & Co., 1927. 
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terminator was about 45” on each night. Other diaphragms were 
inserted along the beam to reduce the amount of scattered moon- 
light. A second observer watched the occultations in the finder, anda 
third operated the oscilloscope and camera. The writer is indebted 
to Dr. James Robertson, director of the Nautical Almanac Office, for 
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Fic. 2.—Observed variation in the light of 8 Capricorni and v Aquarii as they dis 
appeared behind the dark limb of the moon. These are 2-inch sections of a spiral trace 
36 inches in circumference. Time increases from right to left and from the bottom to 


the top. 


furnishing calculations of the times of the occultations from more 
accurate data than those given in the Ephemeris. The predictions 
were accurate to about 2 seconds. The camera was started 20 sec- 
onds before the predicted time, to provide a wide margin of safety. 

The observed curves of light-variation are shown in Figure 2. 
These represent 2-inch sections of the entire 36-inch film. The lower 
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traces in each film represent the time before occultation; the upper 
traces, the time after occultation. Any difference in stability should 
represent the influence of seeing. It is seen that the difference is 
small. This result may be ascribed to the averaging effect of the 
100-inch mirror on the atmospheric waves. 

Although spread out too much in the time co-ordinate, the traces 
immediately suggest diffraction at a straight edge. This possibility 
had not been sufficiently considered by the writer prior to the time of 
observation, although it had been mentioned by Eddington® in 1909 
immediately following MacMahon’s proposal. Very recently J. D. 
Williams’ has published an analysis of the diffraction problem and 
suggested that the small changes in the pattern caused by the finite 
diameter of the source could be used as a gauge of the angular di- 
ameter of stars. The writer’s computation of the expected diffrac- 
tion effects was well along before he learned of Williams’ work. The 
conclusions are the same, and hence the details will not be repeated. 

Williams considers only the case of monochromatic light, and for 
comparison with experiment an extension to the case of a source 
emitting a continuous spectrum is necessary. Figure 3 shows the 
effective spectral response-curve of the photocell obtained by multi- 
plying together the spectral response-curve of the Cs-O-Ag cathode 
for an equal-energy source, the energy-curve for a black body at 
6000" K, and the absorption in the earth’s atmosphere. The altitude 
of 8 Capricorni at the time of occultation was 20°, corresponding to 
an air mass 2.92 times the vertical mass. Only the Go component of 
8 Capricorni is bright enough over most of the wave-length range to 
be effective; the Ao component is therefore neglected. The analysis 
will be carried out in detail only for 8 Capricorni on account of the 
poor signal-to-noise ratio for y Aquarii. The results are qualitatively 
the same for the two stars. 

To reduce the labor of computing the diffraction pattern for many 
wave lengths, the spectral response-curve was divided into five parts 
of approximately equal width. The area and centroid of each part 
was found. The heavy vertical lines in Figure 3 represent the five 

6 M.N., 69, 178, 1909. 


7 Ap. J., 89, 4607, 1939. The writer is indebted to Mr. Williams for the opportunity 
of examining his results prior to publication. 
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monochromatic sources that were computed to represent the whole 
curve. This rather coarse lumping is not as bad as it looks, since the 
parameter which determines the scale of the diffraction pattern is 
proportional to the square root of the wave length. 

Five diffraction-curves were computed by reading off points on a 
standard curve plotted by the Cornu-spiral method.* The five curves 
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Fic. 3.—Effective spectral response-curve of a cesium-oxide photocell, allowing for 


the energy-curve of the star and the absorption in the earth’s atmosphere. The curve 
is lumped into five monochromatic sources, indicated by the heavy vertical lines. 


were then added, and the integrated result is shown in Figure 4. 
This, of course, applies to a point source. To take account of the 
finite diameter of the star, it is sufficient as a first approximation to 
resolve it into an equivalent double with a separation equal to the 
distance between the centers of intensity of the two halves of the 
star.’ For a uniform disk this is 41 per cent of the diameter; darken- 
ing at the limb equal to that of the sun reduces it to 35 per cent. 
For 8 Capricorni the calculated equivalent double star has a separa- 
tion of the order of o%001, which, viewed from the moon’s limb, is a 
distance of 180 cm on the earth. The intensity would then be repre- 


8 See, e.g., Jenkins and White, Fundamentals of Physical Optics, p. 193, McGraw- 
Hill, 1937. 
9 Russell, Dugan, and Stewart, op. cit., p. 745. 
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sented by the sum of twe identical point-source curves, one for each 
component of the hypothetical double, but with one curve displaced 
laterally 180 cm with respect to the other. The resultant is indicated 
by the dotted lines. It has been moved horizontally into the best fit 
with the original curve, because, since the time is not known abso- 
lutely, only the change in the shape of the curve will be observable. 
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Fic. 4.—Comparison of computed and observed diffraction pattern. The telescope 


scans the curve from right to left. 


It is seen that the distortion is not large enough to be observable 
with certainty; if the star had been twice as large, it would have 
produced an effect above the threshold of observation. When the 
separation of the two hypothetical components is small, the de- 
parture from the original curve at the most sensitive points, the 
points of maximum curvature, is proportional to the square of the 
separation. 

The computed noise level indicated an RMS error of about 4 per 
cent of the whole intensity of the star for a single instantaneous 
observation, which checked roughly with the observed values. The 
observed points in the figure have a smaller error than this because 
they were taken from a smooth curve drawn through the original 
record in the same way that a smooth curve is drawn through the 
plate-grain irregularities on a microphotometer trace. The points 
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were fitted by sliding and stretching in a horizontal direction 
only. The observed speed of the fringes past the telescope was 613 
meters/sec, as compared with a computed speed of 541 meters /sec. 
The calculated angle between the moon’s limb and the path of the star 
was 42°9. The observed speed would correspond to an angle of 
50° 5, or a slope of 7°6 on the lunar surface. Such a slope seems not 
unreasonable. The occultation took place at a relatively smooth 
part of the moon’s limb. In the mountainous regions a steeper slope 
might have been encountered. Although it might be anticipated that 
the assumption that the moon’s limb is a straight edge would turn 
out to be too optimistic, the observations on the two stars did not 
show any contradiction of this hypothesis. 

The diffraction-curve is not scanned by an infinitely narrow slit, 
and some distortion might arise from the use of a large telescope. 
At any particular instant the 1oo-inch telescope would receive light 
from a portion of the diffraction-curve 250 cm wide. In order to 
estimate the averaging effect, it is sufficient to replace the telescope 
by an “equivalent double,” as was done for the star disk. The result 
is a separation of about 105 cm on the earth, as compared with 180 
cm for 6 Capricorni; the effect is therefore negligible. With the 200- 
inch telescope, however, the distortion might be troublesome. 

If a potassium hydride photocell as free from microphonic sensi- 
tivity and with as low time-lag as the type of gas-filled cesium oxide 
cell (RCA 921) actually used can be obtained, it would operate with 
a much narrower range of wave lengths. The peaks of the diffrac- 
tion-curve would then be sharper and more sensitive to distortion 
by sources of small equivalent separation. Against this must be set 
the poorer response of potassium-hydride cells to red stars, which 
are, of course, the only favorable cases. Some improvement in 
signal-to-noise ratio could be realized by using a larger value of the 
input resistor. An inspection of Figure 4 shows that the highest fre- 
quency in the observed diffraction pattern was not over 100 cycles 
per second. To record this with negligible phase distortion in the 
amplifier, it would be necessary to have the critical frequency f, be 
of the order of 300 cycles, but not as high as 800 cycles, which was the 
value in the amplifier actually used. A gain of a factor of 2 or 3 is 
possible here. Beyond that the prospect for improvement in ob- 
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servational accuracy seems not very hopeful. It may be remarked 
that the calculation of the signal-to-noise ratio rests on such general 
considerations that it applies to electrometers or to any device 
where the current is sent through a resistance to obtain a voltage. 
There is theoretically a possibility of some improvement by the use 
of secondary emission multipliers, but a trial of two types failed to 
reveal anything as good as the conventional methods of amplifica- 
tion for this particular problem. 

The moon, in all positions of its orbit, occults about a dozen stars 
whose angular diameters, calculated from their magnitudes and 
spectra, are larger’ than o’005. For these stars there appears to bea 
definite hope of obtaining an experimental measure of the diameter 
by observing the diffraction pattern at the time of occultation. 


The writer is indebted to Dr. Joel Stebbins for an accurate calcu- 
lation of the speed of the shadow past the telescope, as well as for 
helpful advice and encouragement; to Mr. A. H. Joy for giving time 
on the 1oo-inch telescope; and to Messrs. D. O. Hendrix and Glenn 
Moore for assistance with the manipulations in the dome. This work 
was supported in part by the Wisconsin Alumni Research Founda- 
tion and the Observatory Council of the California Institute of 
Technology. 

CARNEGIE INSTITUTION OF WASHINGTON 
Mount WILSON OBSERVATORY, 


WASHBURN OBSERVATORY 
UNIVERSITY OF WISCONSIN 
February 1939 


10 For a giant star to have an angular diameter greater than 0”0o05 it must be brighter 
than 1™65 for spectral class Go, brighter than 3™5o0 for class Ko, and brighter than 6™o 
for class Mo. These are visual magnitudes and are calculated from the equation in 
Russell, Dugan, and Stewart, op. cit., p. 612. 








SECOND REPORT ON THE EXPANSION OF 
THE CRAB NEBULA* 


JOHN C. DUNCAN 


ABSTRACT 


The expansion of the Crab nebula, revealed in 1921 by a comparison of the positions 
of twelve nebulous points on photographs made with the 60-inch telescope with an 
interval of 11.5 years, is now studied by observing twenty points on photographs made 
with the 60-inch with an interval of 29.1 years. On the assumption that the observed 
nebulous points started simultaneously and that each has moved with uniform velocity, 
the point from which the expansion takes place is located 20” following and 2” south 
of the northern component of the central double star, and the date of the outburst is 
indicated as A.D. 1172. The observed angular motions, together with the radial veloci- 
ties published by Mayall, indicate for the nebula a distance of about 4200 light-years 
and a major axis of about 6. There is evidence that the nebula has a proper motion of 
0037 per year, corresponding to a velocity across the line of sight of 230 km/sec. 

The south component of the central double star has a proper motion of o%o19 per 
year, while the north component is sensibly stationary. A twelfth-magnitude star 
situated 5/8 north-preceding the nebula has the conspicuous proper motion of 07240 
per year. 





Outward motions of twelve luminous condensations situated near 
the periphery of the Crab nebula (NGC 1952 = M1 Tauri) were 
reported in 1921.' These motions were determined by comparing 
with each other two plates made at the Newtonian focus of the 60- 
inch reflector, one by Ritchey in 1909 and the other by me 11.5 years 
later. Recently I have made similar photographs of this unique ob- 
ject, which, together with the earlier plates, show the expansion 
more clearly. The best of these was made on November 24, 1938, 
with an exposure of 1 hour on an Imperial Eclipse plate; and the re- 
sults now reported were obtained from a comparison of this plate 
with Ritchey’s, the interval being 29.1 years. The measurements 
were made on the blink comparator constructed at this observatory, 
which is superior to the instrument formerly used; and this circum- 
stance, the longer interval, and the superiority of the photograph of 
1938 over that of 1921 combine to make the present results much the 
more trustworthy. 

The co-ordinates x and y refer to rectangular axes intersecting at 

* Contributions from the Mount Wilson Observatory, Carnegie Institution of Washing- 
ton, No. 609. 


1 J.C. Duncan, Mt. W. Comm., No. 76; Proc. Nat. Acad. Sct., 7,179, 1921. 
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PLATE XXIV 
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EXPANSION OF THE CRAB NEBULA 
The vectors indicate the motion of twenty nebulous points in 500 years. Projected backward they suggest an 


outburst occurring about 800 years ago at the point marked by the white spot. Comparison stars are encircled 
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the north component of the sixteenth-magnitude double star near 
the center of the nebula and approximately coinciding with a paral- 
lel of declination and an hour circle. The plate constants were de- 
termined from measurements on twelve stars of the seventeenth and 
eighteenth magnitude, not chosen, as before, near the measured 
nebulous points but situated three in each quadrant and far enough 
from the origin to be free from nebulosity. The probable error of a 
measurement of a comparison star, including any relative proper 
motion, is +o%070 in x and +0%066 in y. The components of the 


TABLE 1 


ANNUAL MOTIONS WITHIN THE CRAB NEBULA 
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central double star were also measured. The north component be- 
trayed no proper motion as great as the probable error, but the south 
component moves in the direction of diminishing right ascension at a 
rate of o’o1g per year. A twelfth-magnitude star, No. 177 on Chart 
402 of the Paris Carte du ciel, situated 4‘0 preceding and 4'7 north 
of the origin of co-ordinates in 1938 (inclosed in a square on Pl. XXIV) 
has a proper motion Ax = —o%058, Ay = —o%229, per year. No 
other conspicuous stellar proper motions were found on the plates. 

Twenty nebulous points were measured, with the results shown 
in Table 1, where the displacements are given in seconds of arc per 
year, and in Plate XXIV, where the vectors represent the motions 
which, at the observed rate, would take place in 500 years. On 
Plate XXIV the comparison stars are encircled. 

On the assumptions that all particles of the expanding nebula 
started from the same small region at the same time and that each 
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has since moved with uniform velocity (though not all with the same 
velocity) relative to the background of stars, the apparent motion of 
each measured nebulous point should conform to an equation of the 
type 

X— % = tds, 


where x is its co-ordinate in 1938 in seconds of arc, x, the correspond- 
ing co-ordinate of the point of outburst, ¢ the time in years since the 
explosion, and Ax the displacement in seconds per year revealed by 
the measurements. From solutions of the twenty equations in x 
and the twenty in y there result the following values: 


%o = +2073 + 272, t = 740 + I09 years, 


Yo = — 158+ 2%0 t= 793 + 23 years. 


II 


The point of explosion thus determined is indicated on Plate XXIV 
by a white dot having a radius approximately equal to its probable 
error. The present center of the nebula is obviously on the oppo- 
site side of the origin. 

The region of the present center is the bright portion in which 
Lampland made the original discovery’ that the nebula was chang- 
ing; but ordinary blue-sensitive plates reveal here no measurably 
definite details and so, while they plainly exhibit Lampland’s 
changes, they leave it‘uncertain whether these changes are due to 
motion or to alterations in the relative brightness. Baade has re- 
cently photographed the Crab nebula at the 1oo-inch telescope, 
using red-sensitive plates, which record beautifully definite details in 
the structure and which will be invaluable for future studies of the 
expansion. On his photographs there appears a bright nucleus about 
8” preceding and 2” north of the origin of co-ordinates. If this is 
the true center from which the nebula originated, it has, according 
to the present. investigation, a proper motion Ax = —o0‘037, 
Ay = +0.005, per year. It is tempting to imagine that this red 
nucleus, although unobtrusive on the photographs, emits the radia- 
tion necessary to keep the Crab nebula luminous but mostly of so 
small a wave length that it is stopped by the earth’s atmosphere. 


2 Pub. A.S.P., 33, 79, 1921. 
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The mean of the two values of ¢t, 766 years, indicates—with con- 
siderable uncertainty—the year A.D. 1172 as the date of the outburst. 

Mayall’s estimate’ of the distance of the Crab nebula, based on 
his value 1300 km/sec of the maximum radial velocity of expansion, 
together with my results of 1921, seems to need but little correction. 
On the assumption that this velocity corresponds to a point 2/7 
from the radiant, where the image of the nebula becomes faint and 
where, according to the present study, the angular displacement is 
o” 21 per year, the distance of the nebula appears to be 1300 parsecs, 
or 4200 light-years. At this distance the proper motion of the nu- 
cleus, found above, denotes a velocity across the line of sight of about 
230 km/sec, and the greatest diameter of the nebula comes out about 
6 light-years. 

CARNEGIE INSTITUTION OF WASHINGTON 


Mount WILSON OBSERVATORY 
February 1939 


3 Tbid., 49, 103, 1937. 








PROMINENCES AND THE SUNSPOT CYCLE 
VINICIO BAROCAS 


ABSTRACT 


The investigations by Lockyer, Evershed, and Bocchino on the relation between 
prominences and spots are extended to all the prominence observations from 1880 to 
1937. A diagram is traced of all the six cycles under consideration, and the behavior 
of high- and low-latitude prominences and their relationship to the sunspot cycle is 
examined. 

The so-called ‘‘metallic prominences”’ are present only in the low-latitude zone. 

The problem concerning the distribution of prominences on the 
solar disk has interested many astronomers. G. Bocchino in a recent 
publication" has continued the works of Lockyer and Evershed and 
has extended the search from 1880 to 1931, about five cycles. She 
has studied the migration of prominences during the sunspot cycle, 
bringing out new characteristics of the latitude and of the frequency 
of the prominences. 

The material already available from the foregoing five cycles ex- 
tended up to 1937 now can be discussed in order to contribute to the 
knowledge of the two zones of prominences (high- and low-latitude) 
and their behavior in connection with the eleven-year cycle. 

Evershed? has already shown that prominences can be divided into 
two different types: low- and high-latitude prominences. The diffi- 
culty, however, was to determine the latitude which separates them. 
In the Evershed diagrams there is a large zone between them, but 
it is so wide that it is impossible to determine a value for the limit- 
ing latitude. Usually, the observers who have worked on this prob- 
lem have taken the mean width of the above-mentioned zone. We 
have chosen the value of + 40° as the separation between high and 
low prominences. Then we will call “high-latitude” prominences 
those which have a latitude greater than 40° and “low-latitude”’ 
prominences those under 4o’. 

We can trace a diagram, as was done by Bocchino,? for the years 

t Pub. Arcetri, No. 51, p. 7, 1933. 

2 Mem. Kodaikanal Obs., 1, ii, 1917. 

3 See also G. Abetti, The Sun, p. 156, London, 1938. 


486 


























PROMINENCES AND THE SUNSPOT CYCLE 487 


1880-1937, taking the years as abscissas and the latitudes of the 
greatest activity as ordinates. These latitudes are determined by 
bisecting the diagram of mean areas at 40° and taking the weighted 
mean of each part. 

The values of prominence areas are taken from Bocchino’s publi- 
cation and the Arcetri’s observations.‘ 

For the low-latitude prominences we write 


We 


SMA TSM B+... $375 % 2 
26 ES Pi 38 


where 2.5, etc., is the mean latitude, and A, B, etc., are the mean 
areas of the prominences; these apply also to the high-latitude 
prominences. 

Table 1 gives the mean values for both the northern and the 
southern hemisphere. The value in the column “Spots” is taken 
for each year from the Monthly Notices. From this table the dia- 
gram in Figure 1 has been traced; it confirms Lockyer’s and Ever- 
shed’s deductions. 

Let us now analyze the curve of low-latitude prominences in the 
diagram. It follows the spot development very closely, although it 
is limited to a zone between + 18° and + 28°, while the spots are 
found, as is well known, in a zone between + 7° and + 23°. 

Prominences appear at a distance of 10° north and south from the 
spots and follow them, but with a smaller range of latitude. At the 
minimum period of the sunspot cycle the latitude of the low-latitude 
and spots is a maximum. The minimum latitude for 


prominences 
and spots is reached one or two years before sunspot 


prominences 
minimum. Yet this law is not a rigid one; it is sufficient to observe 
that, when spots reach their minimum of latitude, the prominences 
of low latitude, while continuing their motion, are moving toward 
higher latitudes after attaining the minimum just referred to. Dur- 
ing the period when the spots reach their lowest latitudes, we have 
the greatest distance in latitude between spots and low-latitude 


prominences. 
The determination of the correlation coefficient will confirm what 


4 Pub. Arcetri, Nos. 52-56, 1934-1938. 
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has been said about the analogy which exists between the curves of 
the low-latitude prominences and of the spots. Mineur’s method has 
been followed in this investigation, and from Figure 2 it is found that 
the correlation coefficient is about 1 (cos ¢ = 0.99). In Figure 2 the 
two straight lines start from above ro” on the abscissa for, as is clearly 


TABLE 1 
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1880. 25-5 54°1 19°7 |} 1909 20°7 49°0 9°7 
POOL 5:5 23.8 57.9 18.2 | IgIoO ; 19.0 52.7 9.9 
1882.. 223 64.8 17.6 | TOLL: . : 23.2 52.6 6.9 
1883. . 22-3 54.38 12.4 IgI2 22.0 49.8 14.1 
1884. 19.8 67 FE.2 fl FOUR... 4. 28.4 49.4 224 
1885. . 20.4 50.4 II.5 |} 1914 20.8 49.8 22.7 
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1887. . 23..5 50.8 S:5 ] 1916 23.9 50.7 16.0 
1888. . 23.0 Se * eee: IQI7 PD eae O1.1 14.2 
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1897.. 18.7 52.9 8.0 1926 23.4 62.1 18.6 
1898. . 20.4 50.6 10.3 1927 25.3 61.8 11.6 
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IQOl.. 20.2 53-4 12.4 1930 20.9 50.0 9.8 
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seen by the dispersion of the points between 6° and 10°, if their mean 
is taken, there will be two perpendicular lines, i.e., a correlation 
coefficient of o, the angle o being go” (cos ¢ = 0). 

Let us now see what the meaning of this result is. A correlation 
of form undoubtedly exists between the two diagrams under dis- 
cussion; but it is almost zero for the smallest value of latitude that 
can be deduced from what has been said, and particularly by con- 
sidering the years 1888, 1911, 1922, and 1932 in Figure 1. 
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Turning now to the diagram of the high-latitude prominence, we 
see that it is more difficult to find an analogy between this curve and 
that of the spots. Prominences drift toward higher latitudes and 
reach them about the time of sunspot maximum. Then, while the 
high-latitude prominences are, in a given year, at their highest lati- 
tudes, low-latitude prominences in the same year reach their lowest 
latitudes. The low-latitude prominences attain the highest latitudes 
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when the high-latitude prominences are at a minimum. It is rather 
difficult in this case to find the correlation coefficient, since the two 
diagrams are not in phase. A remarkable result is obtained by dis- 
placing the diagrams of the high-latitude prominence by three years 
with reference to that of the spots. This gives an almost perfect 
coincidence between the maxima, but the minima are unaffected by 
the displacement: the minima of the prominences precede those of 
the spots. The dispersion of the points in Figure 3 shows that there 
is no correlation in form between the two curves under considera- 
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tion. The two lines are almost perpendicular in low latitudes, and 
they form a small angle in high latitudes; in the first part there is no 
correlation; in the second one should naturally exist, but it is a 
forced correlation, for the curve has been displaced in order to ob- 
tain the coincidence of maxima. It can therefore be stated that a 
correlation does not exist between the curve of the high-latitude 
prominences and that of the spots, or at least only as far as it con- 
cerns the duration of the cycle. 

In Figure 1 we see that the migration of high-latitude promi- 
nences toward the poles is not connected with the phenomena of the 
low-latitude prominences. Yet they reach a maximum latitude when 
the spots and the low-latitude prominences drift toward the equator, 
that is to say, some years after the spots and the low-latitude promi- 
nences have reached their maximum latitude. This fact occurs, as 
is well known, about the period of sunspot maximum. During the 
maximum of the cycle we have, then, a greatly disturbed zone of low 
latitude at a mean latitude of about + 20°, while at the high- 
latitude zone the prominences are at their maximum latitude of 
about 60° and are already decreasing in activity. 

During the period of sunspot minimum the low-latitude promi- 
nences also appear at their highest latitudes; this fact may be at- 
tributed to the presence of two zones of activity. By a careful study 
of these prominences and of dark flocculi a distinction can be made 
between the prominences of the old and the new cycles; at present, 
except for the latitude, there is no criterion by which we can make 
this distinction. In the case of the spots, the sign of the magnetic 
field provides a definite distinction. 

The present investigation may also explain why at the minimum 
period there are two zones of activity at low latitudes and only one 
at high latitudes. 

To give a general view of the phenomenon, another diagram (Fig. 
4) has been drawn in which the abscissas represent years and the 
ordinates represent the mean areas expressed in prominence units. 

The curves representing the motion of the prominences (low- and 
high-latitude) are the means for the northern and southern hemi- 
spheres for the six cycles from 1880 to 1937.5 The ordinates of the 


5 Values taken from various issues of Pub. Arcetri. 
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curve for the spots are Wolf numbers taken from the Zurich publi- 
cations. In Table 2, their values are given. The prominences are 
divided, as before, into low- and high-latitude prominences. 
Analyzing Figure 4, it will be noted that in this diagram the low- 
latitude prominences follow the spots, i.e., they follow the spots not 
only in the latitude distribution but also in frequency. High-latitude 
prominences, on the contrary, behave in a different manner: they 
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reach a maximum about two years before the maximum of the sun- 
spot cycle. 

To give a clear-idea of the phenomenon, the diagram in Figure 5 
shows the curves representing the means of the six cycles given in 
Figure 4 for spots and for high- and low-latitude prominences. In 
Figure 5 the abscissas represent the years of the cycle and the ordi- 
nates represent the mean value given in Table 3, expressed in promi- 


nence units. 


years 








PROMINENCES AND THE SUNSPOT CYCLE 493 


The first thing to be observed in Figure 5 is that the curve for the 
spots reaches its maximum in the fifth year. 

High-latitude prominences reach a maximum two years before the 
spots and then decrease to a minimum five years before the spots; 
they have, therefore, no connection with the spots. But the cycle 




















|| 
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1880 20%..2 262.9 42.9 | 190g 3 ye 82.5 43.9 
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of the high-latitude prominences also lasts eleven years, although 
it is not in phase with the spot cycle. Its maximum appears in the 
eighth year, if we take the minimum of the curve as the beginning 
of the cycle, as in the case of the spots. 

More interesting is the curve of low-latitude prominences, which 
behave differently. They, too, reach a maximum two years before 
the spots, but afterward they decrease and reach a secondary mini- 
mum one year before sunspot maximum. Then they increase and 
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attain another maximum, higher than the preceding, at the same 
time that the spots do. The minimum occurs in the same year as 
that of the spots. The curve of the low-latitude prominences follows 
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very nearly that of the spots with respect to its maximum, its mini- 
mum, and its amplitude, except for its behavior in the third year 


before sunspot maximum. 
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From Figure 5 we see that during the period of sunspot minimum 
the low-latitude prominences also have a minimum while the high- 
latitude prominences are already advanced in their cycle. At the 
beginning of the spot cycle the three curves increase rapidly; and 
the high- and low-latitude prominences reach a maximum in the 
same year, exactly two years before the spot-cycle maximum. The 
high-latitude prominences then begin to decrease and are followed 
by the low-latitude prominences. One year before the spot-cycle 
maximum the activity of the spots, which has reached a very high 
value, begins to make itself felt in the low-latitude prominences, and 
the latter suddenly increase and reach a maximum at the same time 
as the spots; then they follow the behavior of the spots, up to their 
respective minima. 

The high-latitude prominences, after the secondary minimum of 
the low-latitude prominences, slowly decrease; but from the time of 
maximum activity they hasten toward their minimum. 

It has been shown by G. Abetti® that the shape of the corona de- 
pends upon the development and the frequency of the prominences 
in high latitudes. 

Having distinguished between high- and low-latitude promi- 
nences, it is interesting to investigate the difference in their consti- 
tution. Several observers have classed them as metallic and hydrog- 
enous prominences according to their spectra, because in the former 
many metallic lines are present, while in the latter only hydrogen 
lines appear. 

The Kodaikanal Observatory publishes twice a year a list of the 
metallic prominences, with their positions, their height, and the 
origin of their lines. The metallic prominences usually differ from 
the others by the presence of b,, b., and b, of Mg; D,, and D, of Na; 
and a group of lines of Fe 1, Fe 1, He, and Cr. I have examined these 
lists from 1912 to the first half of 1938 and have noticed that the 
prominences present more or less the same lines, i.e., those given in 
Table 4. Several prominences show numerous lines of other ele- 
ments; but if we take their height, we find that they are so low that 
they may be regarded as more like the chromosphere itself than as 


® Ibid., No. 56, p. 53, 1938, and Mem. Soc. Astr. It., 9, 161, 1938. 
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prominences, and this probably explains the presence of lines that 
are not usually found in prominences. We define ‘‘prominences”’ as 
all eruptions that attain a height of 20’’ (about 14,000 km—the 


height in the chromosphere of Ca 11) or more. 
The latitude at which metallic prominences appear is what we 


TABLE 4* 


Chromosphere | Elements Id If Height | E.P 
4923.9 ; Fe Il : } 30 2000 2.88 
5015.6 He 2 2500 20.52 
5018.4 Fe il 4 25 2000 2.88 
5167.3 Meg IS 18 1500 2.40 
5167.5 Fe 5 18 1500 1.48 
5168.9 Fe 7 25 1500 > 88 
5171.6 | Fe 7 25 1500 2.88 
5172.6 Mg 20 30 2000 | 2.70 
phe Meg 30 40 | 2500 2.70 
5197.7 Fe Il 2 12 500 2 22 
5204.8 Cr-Fe 8 7 | 500 0.04 
5206.2 Cr 5 8 600 0.94 
5208.4 | Cr 5 12 600 0.04 
5227.0 | Fe-Cr | 8 8 500 b.55 
5234.6 | Fell 2 15 | 500 | 3.21 
5269.5 Fe 8 15 600 0.86 
5270.6 |  Fe-Ca 7 7 500 1.60 
5275.9 | Fell 6 | 20 500 | 3.19 
5284.1 | Fell I 8 400 2.88 
5316.6 Fe Il 6 30 850 3.14 
5324.2 Fe 7 i 400 | 3.20 
5302.8 Fe ll ‘< 15 500 | 3.19 
= 382 3 Fe 6 6 400 | 

5889.9 Va 30 25 1500 | 0.00 
5895.9 Na 20 20 1500 0.00 
6678.1 He-Fe 5 20 2200 21.13 
7005.2 | He 6 1000 20.87 





se Pe intensity on the disk; J, = intensity in the flash. 


have called “low latitude,” i.e., less than 40°. Table 5 gives a list 
of the metallic prominences in high latitudes observed in the period 
1912-1938. We find that during a period of 26.5 years, out of 835 
metallic prominences observed on the sun, only 28, or 3 per cent, 
were in high latitudes; and all showed the same lines, namely, b,, b., 
b,, D,, and D,. 

When we take all prominences, irrespective of whether their 
spectra show metallic or hydrogen lines, the percentage of high- 


latitude prominences is 41. 











Date 


IQI2 
Mar. 19 
Sept. 31 

1913 
Mar. 13 
Mar. 26 
Mar. 26 

1914 
April 6 
Get... 11 
Nov. 4 
Nov. 26 

IQL5 
Jan. 15 
May 20 

1910 


Jan. 


IQIgQ 
Jan. 28 
Feb. 23 
Mar. 2 
Mar. 2 

1920 
Jan. 
Jan. 
Feb. 
Nov. 30 
Dec. 3 
Dec. 12 

1921 

1922 

1923 

1924 

1925 

1920 
Feb. 16 
Feb. 17 

1927 
Mar. 24 

1928 

19209 

1930 
Dec. 26 

1931 

1932 

1933 

1934 

1935 

1930 

1937 

1938 


I2 
17 
5 
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TABLE 5* 


Lat 
52.015; 
78.5 N. 
46.5 S. 
44.5 8. 
4. & S. 
45::0°9. 
49.5 S. 
50.0 N. 
52.0N. 
53.6. 
Ca ts: 
42.058. 
42.0.9: 
41.058. 
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Na 
Na 


Na 
Na 
Na 


Na 
Na 
Na 
Na 


Na- 
Na- 
Na- 


Na 
Na 
Na 


Mg-Fe 


Mg-Fe 
M g-Fe 
M g-Fe 
-M g—Fe 


Elements 


Mg-Fe 
Mg-Fe 


Mg-Fe 
M g-Fe-He 
M g-Fe-He 


Fe 
Fe Il 


Mg-Fe 


Mg-Fe 
Mg-Fe 


M g-Fe-Fe il 
Mg-Fe 
Mg-Fe 


Fe il 
Fe it 


Fe il 
Fe Il 


Mg-Fe 
M g-Fe 
Mg-Fe 
Mg-Fe 
Mg-Cr 
M g—Fe 


Fe il 


Fe 11-Cr 
Fe 11-Cr 


M g—Fe-He-Fe i 
-M g—Fe 


M g—Fe-He-Fe il 





_ ,* The second column gives the number of the metallic prominences for the year. Values for 1938 are 
limited to the first half-year. 
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The difference between eruptive and quiescent prominences is 
perhaps only a difference in intensity; Pettit has shown that the 
spectra of quiescent prominences contain lines with an intensity 
greater than 30 on Mitchell’s scale, while the intensity of the faint 
lines of the metallic prominences is 15 on the same scale. 

The Kodaikanal observations in the region between 4923.9 A and 
7065 A show eleven lines in the metallic prominences fainter than 
15 in Mitchell’s scale (see Table 4); these lines include those of C7, 
Fe, and Ca. The three lines of ionized barium which have intensities 
of 25, 20, and 20 do not appear. 

What has been said about the metallic prominences observed be- 
tween 1912 and the present time agrees with Pettit’s examination’? 
of the plate of the eclipse of June 8, 1918, with the exception of some 
lines fainter than 15 observed at Kodaikanal. In Table 4, Ja is the 
intensity on the disk and /; in the flash, according to Mitchell. 

Whether the differences in the spectra of both types of promi- 
nences depend only upon the intensities of the lines could only be 
decided after a systematic study of the spectra of prominences photo- 
graphed by Lyot’s method. The eruptions are certainly more vio- 
lent at low latitudes, where the excitation of the elements in the 
chromosphere and prominences must be greater than in the high- 
latitude zones. 

R. OBSERVATORY ARCETRI-FIRENZE 


November 1928 


7 Mt. Wilson Contr., No. 451, p. 2; Ap. J., 76, 9, 1932. 








THE POINT-SOURCE MODEL WITH 
CONSTANT OPACITY 
PHILIP C. KEENAN 


ABSTRACT 

The differential equations for the point-source model, in the form derived by 
Chandrasekhar and von Neumann, have been integrated numerically for eight values 
Mw ranging from 19.3 to 31,500 in units of the sun’s mass. All the solutions are char- 
acterized by negative density gradients in the central region, since p becomes zero at the 
center. The radius at which the density is a maximum shifts outward with increasing 
mass so that é*, the fraction of the radius containing 90 per cent of the total mass, in- 
creases steadily from 0.65 to 0.98. For large masses the central regions are practically 
isothermal. 

The mass-luminosity curve becomes less steep with increasing mass, approaching 
the line L~M. 

The physical interest of the point-source model is increased by the agreement found 
between the model and the massive Trumpler stars, both in the values of &* for corre- 
sponding masses and in position in the mass-luminosity diagram. 

The stellar model with the sources of energy concentrated effec- 
tively at the center forms an important limiting case and was con- 
sidered first by Eddington,’ who obtained a numerical solution for 
one set of values for M, L, R, and uw. Further discussions have been 
given by Cowling,? Milne,* Stroémgren,* and Biermann,’ but the 
most complete development of the general properties of such models 
is due to Chandrasekhar and von Neumann.° 

Among the point-source models, those with constant opacity (x) 
provide a set of limiting configurations in so far as at very high 
temperature the main source of opacity arises from electron scatter- 
ing (for electron scattering x = const.). The complete discussion of 
such models requires considerable care, for, if we start with the as- 
sumption of radiative equilibrium at the surface and continue in- 
ward, a stage will be reached where the radiative gradient will exceed 
the adiabatic gradient. Here the radiative gradient will become un- 
stable, and it is generally assumed that beyond this point the 

1M.N., 85, 408, 1925; The Internal Constitution of the Stars, p. 124, 1926. 

7M.N., gt, 92, 1930. “ZF. AD, 2s 3455. 193%: 

3M.N., 91, 4, 1930. 5 Zs. f. Ap., 3, 116, 1931; tbid., 4, 61, 1932. 

®S. Chandrasekhar, An Introduction to the Study of Stellar Structure, chap. ix, p. 332, 


1939. 
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adiabatic temperature gradient will be set up. Although a rigorous 
discussion has been given only for the case where the radiation pres- 
sure as a factor in hydrostatic equilibrium is negligible, there can be 
no doubt that a continuation of the equations beyond the region of 
known stability requires caution. In spite of this, the establishment 
of the (L, M, x, u) — relation in detail for the complete point-source 
model is of particular value in that it should provide a ‘minimal 
case.” 

In view of these considerations, numerical integrations have been 
carried out for the compiete point-source model with « constant, in 
which the mass of the star has been exactly used up when the center 
has been reached. No attempt has been made to deal with models 
with point mass at the center. 

The equations of equilibrium, 


d(py + pr) _ _GM(r) 


dr r? 
and 

dp, KL 

; P eaCeTee reget a 

dr 4mcr? 


are transformed by the substitutions:’ 


yr=atk=-; P, = Iz ; Pe = HZ, 
x 


where the constants are connected by the relations 


ne KL Rk (3 1/4 ts v 
* | 16%¢G pH = : 


4mc? k 4 3 1/3 
= — ~--—— a II, 
= | ono (5) 3 





into the forms 


dx ; 
7. rer 


7 Cf. ibid., for notation and derivation of equations I and IT. 
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and 
d?z 
—_ = — 4 
‘Tr x (II) . 


For the case under consideration the three necessary boundary 
conditions are: 


a) At the surface, where? = 0, Z=0. 
dz 

b) At the surface, where? = o, = o> 'O 
C 

c) At the center, where ¢ = ¢. , Mrao = 0. 


Furthermore, since equations I and II admit of a constant of 
| homology so that the equations are invariant to the transformation 


x, 2,t—>C4x, Cz, Ct 


where C is an arbitrary real number, we can normalize the units of 
« so as to make x = 1 fort = 0. 

It is shown in Chandrasekhar’s discussion® that the solutions satis- 
fying these boundary conditions are those for which the zé-curve 
passes through a maximum and intersects the taxis at ¢, with a 
slope (dz/dt). = —6,. The mass and luminosity relations for these 
models are 


M = 1.117 © 86)/? (6; + 62)u~? 


(GM 5 I 
«oo aT m= 1.471 X 10% « "yu * 63/* . 
K I 2 


L 

The constant 6, is uniquely determined by 6,, and it can be shown 

| that for a sufficiently large positive value of 6,, 6, must be positive. 
| As 6, decreases, 6, also diminishes, and for a critical positive value 
(6,)., 5. becomes equal to zero, corresponding to a model of zero mass. 
The general procedure in carrying out each solution is to start at 
the surface with some value of 6,, presumably greater than the as yet 
unknown (6,),, and to integrate inward until z = o, where /, and 6, are 
then found. At this point r = o, and, hence, x goes to infinity. As 6, 
approaches the critical value, the range of integration becomes great- 
er and the numerical work becomes increasingly difficult, but this is 
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not serious, since we are primarily interested in models of large mass 
(large 6,), for which no solutions have been attempted previously. 

In the actual numerical work it is convenient to introduce new 
variables 7 and ¢ defined® by 


Then I and II become 


dx = §-3/4¢ 1-1 4 (1’) 
dr 

i eo z 
dr? = LA (II’) 


and the initial conditions at the surface of the star are now: 


dé 


Atr=o0; 7-7, f=0, —-=<= 1, 
dr 


Since dx/dr is infinite at the boundary, it is necessary to obtain 


series expansions to start the solutions. Thus, for x and ¢ near rt = o 
we find by repeated substitution and integration: 


ge ge Sy ge 4 4 Set , . .) 


= 6-3 a (51273 © al 332775 4 404, 3527/2 . ee 


135 904, O81 
0/4 (5127/4 133,072 5/2 10, 496 11/4 
+ 6; ii + 10,125" + 1155) Pissed 
6 BGs 
2 


a si ao 256 } 128 2 64 - 
Gag ae mm et et a ee ae i a 


5 


= 5: 3/2/1285 ® she 25»472,7/2 4+ 79/74 .., a 
3 4725 
PO ae Re sees) 


For small masses the series converge slowly, and it was found more 
convenient to extend the first one or two values by Kutta’s modifica- 


8 For the suggestion of this transformation I am indebted to Professor Walter Bartky. 








THE POINT-SOURCE MODEL 503 


tion of Simpson’s rule to the point where ordinary forward integra- 
tion could be begun, than to work with extremely small intervals 
in T.° 

Solutions were made for eight values of 6,, corresponding to a 
range of 19.3 to 31,500 in My’, where M is expressed in units of the 
sun’s mass. From a plot of 6, against the corresponding values of 
5, + 6., it was then possible to estimate that the critical value 
(5,)o, for which 6, = 0, is approximately 8. 

In preparing Table 1, which gives the results for the individual 
models, the solutions have been condensed considerably, particularly 
in the case 6, = 11.51, for which the original integration required 229 
steps. The second column of the table lists the temperature ex- 





TABLE 1 
6; =625 M wp? =31,500 | 61 =107.3 M yp? =19054 
ry 47.225 (£/T4) max =1.019 | tr 1=1.3420 (¢/7 4) max =0.6650 
r/R=1/x T/T¢ P/ Pmax | r/R=1/x T/T¢ P/ Pmax 
©.ggo 0.258 0.050 0.963 0.236 0.047 
985 . 386 .190 946 .352 0.156 
.982 459 311 936 .419 0.256 
979 . 546 498 924 498 0.414 
970 .604 .640 gio 593 o.640 
974 649 752 goo .656 0.796 
Q72 .686 .838 8oI .705 ©.QOol 
.97I 718 gob .884 745 °.g66 
969 747 953 877 780 ©.996 
968 E72 983 $73 796 I .000 
966 795 997 870 SII ©.9g90 
965 .816 .998 863 838 ©.970 
964 .836 g86 848 . 886 0.846 
962 854 .g60 831 928 0.640 
961 872 .920 807 .964 0.3605 
059 888 .870 . 786 O81 ©. 206 
.958 .9O3 809 769 9890 ©O.122 
956 g18 .736 0.730 0.997 0.036 
954 .932 .053 
952 946 558 
.949 .958 454 
940 .Q7I . 340 
O41 .983 .218 
0.932 ©.9904 0.087 




















9 Cf. Levy and Baggott, Numerical Studies in Differential Equations, I, 106, 1934. 
The forward integration was done by their Method III of chap. iv. See also pp. 159-60. 








TABLE 1—Continued 





























=O 


6:1 =73.1 
Th 4=1. 4013 (¢/r' *)max 
7/R=1/x T/T¢ 
0.952 0.220 
929 0.338 
.Q17 ©.401 
go2 0.477 
884 0.508 
.871 0.6025 
861 0.675 
852 0.714 
836 0.770 
828 0.803 
821 0.827 
813 0.849 
805 0.870 
797 0.888 
779 0.923 
755 0.954 
738 0.9090 
714 0.983 
0904 o.ggo 
650 0.990 
0.135 I.000 
61 = 42.6 M yp? = 383 
7 41.537 (¢/r" *) max 
r/R=1/x T/T, 
0.9290 oO. 200 
898 308 
. 880 . 300 
860 435 
.830 517 
$20 S73 
500 O15 
.790 OSI 
oY | . 708 
. 700 754 
-745 -792 
- 730 520 
-715 556 
.699 883 
.682 QOd 
.663 931 
640 953 
608 973 
545 992 
O-Si3: 0.990 








SP9S9O0DDDDD0D0OOHOCOOOOO0O0 0 
roe 
NO 


8oqI 


P/ Pmax 


0.004 
O19 
onus 
314 
490 
633 
741 
524 
930 
Qggo 
9908 
g06 
QO5 





0.90 
$68 
540 
822 
-793 
774 
-759 
747 
.720 
7O9 
094 
673 
654 
635 
.616 
597 


t 
Num s] 
nn oO 


hpAAUMY 
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o+ 
Wm NW 











.9388 
0.999 





M p?=184 


(¢/r' 


0.904 


4) 
max 





0.806 


=I 


P/ Pmax 


C©C0O0O0O0C0O0C0CO0OF COO OOOO 0 0 
+ 


an 


102 


P/ Pmax 


0.029 
0g4 
150 
250 
409 
526 
620 
. 700 
$20 
O13 
g00 
999 
057 
Q40 
503 
764 
.043 
. 508 
302 
. 208 
102 
0.050 
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TABLE 1—Continued 











61 =16.90 M yu? =56.0 b: = 11.51 M p?=10.3 
r 4 = 2.167 (¢/r" ‘) max =1-778 Te “=3.020 (¢/r* ‘) max =3 174 
r/R=1/x T/T¢ p/Pmax r/R=1/x T/T¢ | P/ Pmax 
0.868 0.146 0.018 0.832 0.105 0.010 
814 .218 0.059 806 .¥2 .O17 
786 . 260 0.097 By vo .148 028 
755 . 309 0.161 734 186 .055 
720 307 ©. 260 698 242 .Og! 
607 . 406 0.342 670 245 322 
680 430 O.412 646 278 .174 
660 .402 0.473 025 303 . 219 
644 502 0.576 602 331 .278 
627 530 0.661 575 . 366 . 301 
605 574 0.762 555 394 -431 
582 617 0.859 535 -420 -§02 
562 653 0.926 515 .450 .582 
545 683 0.968 104 480 .664 
529 710 0.992 468 .518 . 766 
514 734 I .000 443 557 856 
. 500 .756 0.994 403 616 .958 
487 «770 0.979 472 662 .997 
474 795 0.953 340 . 700 .990 
449 829 0.880 313 .748 .948 
424 859 0.785 280 793 857 
399 . 886 0.6075 248 .833 -739 
373 gil 0.553 216 872 597 
.344 .934 0.426 177 .QI4 .418 
312 .956 0.206 131 .958 ael 
271 .976 0.164 .104 .978 539 
244 .986 0.100 082 .QgQo 049 
200 995 0.035 0090 995 024 
0.120 0.9995 0.004 0.O5I 0.999 0.005 











pressed as a fraction of the central temperature 7., and the third 
column gives the density in units of the maximum density occurring 
in the given model. These relative units are used because of the 
homology invariance of the models considered. However, the scale 
factor*can be determined for any assigned value of R, and by mul- 
tiplying the quantities in columns 2 and 3 by the appropriate factors 
the physical variables can be expressed in any desired units. For this 
purpose the values of ri/* ~ T, and (£/7'/*)max ~ Pmax are given at 
the top of the table for each solution. 

The density and pressure distributions for six of the models are 
presented graphically in Figure 1. 

We are now able to determine the slope of the mass-luminosity 
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relation. To put it in the usual form, the equations given above for 
M and L are written 


log M = log 1.117 © — 2 log uw + log 6!/7(6, + 6,) 


I.47I X 104 


oi — 2.5 log 63’. 


2.5 log L = 2.5 log 


When log My? is plotted against 2.5 log 63/7, as in Figure 2, we 
obtain the ordinary mass-luminosity diagram except that the curve 





























T q t ‘ q t ' ' 1 
54 
i 19.3 56.0 184 600 se” 31,500 
0.8 ~ 
0.60 IP, Pmax 
0.4Fr . 
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1.0 31,500 
1954 
°o.8F 600 aa 
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56.0 
— 19.3 5 ee 
0.4 
O. 2b 
l l l | 1 l l ] l 
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Fic. 1.—Distribution as density and temperature in the point-source models. 
Abscissae are given as fractions of the total radius. 
is displaced both horizontally and vertically by constants depend- 
ing upon the values assigned to uw and x. For large masses the slope 
of the curve is much less than for small masses, where the point- 
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source model is known to agree with the Standard Model in pre- 
dicting L ~ M**%, approximately. 

The computed values are collected in Table 2, in which the last 
column represents the quantity &* defined by Chandrasekhar’ as 
that fraction of the radius which incloses 90 per cent of the total 
mass. It is thus a measure of the degree of central condensation of 
the matter of the star, and is about 0.43 for main-sequence stars. 


2.5 log 63/2 
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Fic. 2.—Mass-luminosity diagram for point-source models. Ordinates are propor- 


“. 


tional to absolute magnitudes. The broken line below the main curve represents the 


slope corresponding to L~M. 


For our model ~* can be found by substituting ds/d¢ for 6, in the 


mass equation, which then runs 


pee ye dz a 
M; = 1.117 Op 76; (F +6.) 


Dividing by the expression for M and setting 


M; 


aa ‘ ate * 
yw °-9 for & = &*, 


we find 
(3 0.9 (6, + 4:) — 6 
= = o€ I a) ~~ G3. 
dt ee : 


10 Op. cit., chap. vili, p. 292. 
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This expression gives (dz/dt);« and hence (d¢/dr):*, so that &* can 
be obtained by interpolation from the numerical solution for the 
given model. 

For actual stars £* can be found when M, L, R, and uw are known, 
and, since the massive Trumpler stars are known to have a relatively 
low degree of central condensation, it is of interest to note that the 
values determined for them by Chandrasekhar" have the same gen- 
eral run as those computed for point-source models of the same order 
of mass. This agreement is merely suggestive, for, as Chandrasekhar 


























TABLE 2 

br 52 M yp? Lyx X10 4 log My? log 63/2 t* 
II.51 71 19.3 0.328 1.29 0.348 0.652 
10.9 S25 56.0 Lage L.75 1.068 725 
29.0 14.4 184 8.00 2.26 1.74 .805 

2.6 25.4 383 18.85 2.58 270% 852 
54.3 35.6 600 S52 2.78 RE 873 
73-1 52.7 1,020 50.3 3.00 2.58 597 
107 .3 84.9 1,054 111 3.20 2.88 922 
025 5606 31,500 1,970 4.50 4.13 0.979 





| 
| 
| 
| 
| 
| 
| 
| 
| 
| 
| 





has pointed out,’ any model characterized by sufficiently extensive 
negative density gradients will give large values for &*. On the other 
hand, it may be further noted that, if the mass-luminosity relation 
for the point-source model is fitted to the observational curve for 
ordinary stars by extrapolating it to small masses, it represents very 
well the position of the Trumpler stars in the diagram.” 

In considering the possible physical significance of this model, it 
must be emphasized again that the properties derived represent an 
extreme limiting case in the discussion of the mass-luminosity re- 
lation. 

These computations were undertaken at the suggestion of Pro- 
fessor Chandrasekhar, and it is a pleasure to acknowledge helpful 
discussions with him at all stages of the work. 

YERKES OBSERVATORY 
January 1939 
1 [bid., p. 314. 2 M.N., 97, 149, 1936. 


13 For the observed mass-luminosity relation see Kuiper, Ap. J., 88, 489, 1938. 





























ON THE COLOR OF P CYGNI 
FRANCES SHERMAN AND W. W. MORGAN 


ABSTRACT 


From a comparison of color excesses of stars near P Cygni and located at about the 
same distance from the sun, it is found that the color excess of P Cygni is partly due to 
interstellar reddening and partly to conditions in the immediate neighborhood of the 
star itself. 

Attention is called to a systematic progression of determinations of color excesses 
apparently depending on the range of wave lengths observed. 


I. INTRODUCTION 

At the present time the evidence seems conclusive that most, if 
not all, of the reddening observed in distant ordinary B stars is 
interstellar in origin. In the case of the B and A supergiants there 
may be intrinsic differences of a minor amount from dwarfs of the 
same spectral type, but the greater part of the reddening found for 
these objects also seems definitely interstellar in nature. 

The B stars having strong emission lines, however, seem to be- 
have differently. They were found by E. G. Williams’ to be system- 
atically too red for their distances as determined from the inter- 
stellar K line in their spectra. He suggested, as the probable explana- 
tion, that their light is subjected to selective scattering in an exten- 
sive atmosphere surrounding the star; the existence of such extensive 
atmospheres had been predicted earlier by Struve from his theory of 
the origin of bright lines. 

Further evidence on the subject has been found recently from 
observations of y Cassiopeiae. Chalonge and Safir,? and Greaves and 
Martin,’ found that the color temperature of this Be star decreased 
considerably during its recent light and spectral variations. From 
gradients determined by the Greenwich observers the temperature 
was found to be 16,100° in 1926 and 9g,200° in the latter part of 
1937. The first temperature corresponds to that of a star of type B8, 
while the latter is similar to one of type A5-Fo as determined at 
Greenwich. In the case of the changes in y Cassiopeiae there can be 
little doubt that they were due to conditions in the immediate 

t Ap. J., 79, 303, 1934. 7 C.R., 203, 1329, 1936. 3 M.N., 98, 434, 1938. 


509 








510 FRANCES SHERMAN AND W. W. MORGAN 


neighborhood of the star; we thus have an individual case in which 
selective absorption has been shown to exist as an atmospheric 
phenomenon in a Be star. 

It is the purpose of the present paper to attempt to find what part 
of the color excess of P Cygni is due to ordinary selective absorption 
and what part—if any—is peculiar to the neighborhood of the star 
itself. This was done by determining the color excesses of stars near 
P Cygni in the sky and situated at the same distance from the sun. 


II. THE DISTANCE OF P CYGNI 

Spectroscopic determinations of the luminosity of P Cygni are 
of doubtful value, since most of the observed absorption features do 
not originate in the star itself but in the expanding envelope. The 
safest procedure is probably to use the intensities of the interstellar 
K and D lines as criteria of distance. For the D lines the intensities 
and reduction curve by Merrill* were used, and an apparent dis- 
tance modulus of .9.1 was found. This was determined from a 
general distance-intensity relation for all stars, as regional differences 
had been found by Merrill to be fairly small. 

A different procedure was followed in the case of the K line. The 
measured total absorptions determined by Williams® (unreduced to 
an exponential contour’) were plotted against apparent distance 
moduli for stars for which Williams had determined spectroscopic 
absolute magnitudes.’ Stars from galactic longitude 340° to 85 
were included which are within 10° of the galactic equator and 
which have spectral types between Bo and Bs. The Be stars were 
omitted. From the K line an apparent distance modulus of 9.5 was 


found. 
III. THE COMPARISON STARS 


Fifteen stars located within about half of a degree of P Cygni were 
selected to serve as comparisons. For nine of these having spectral 
types between B8'and A5 the Draper types were used, and the as- 
sumption was made that they are main-sequence stars. Each star 

4 Ap. J., 86, 297, 1937; 87, 120, 1938. 5 [bid., '79, 291, 1934. 


6 Irregularities in the observed contours of interstellar K and D lines, as found by 


Beals, make the use of the uncorrected total absorption preferable. 


7 Ap. J., 83, 333, 1936. 
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was assumed to have the absolute magnitude of the mean of the 
main-sequence stars of its type, as estimated from stars in the 
neighborhood of the sun. Spectra of the other objects were obtained 
at the Yerkes Observatory with low-dispersion cameras attached to 
the Bruce spectrograph. 

The luminosities of the B stars were based principally on the in- 
tensities of the hydrogen lines, while the luminosity of the F5 
dwarf was determined from the enhanced lines commonly used for 
stars of this type. The luminosity standards for the early-type stars 
were y Orionis (B2) and x? Orionis (cB1); the calibration of giants, 
intermediates, and dwarfs by Williams was used. The absolute 
magnitude of y Orionis was found to be —3.6 by Kapteyn’® in his 
study of the parallaxes of the stars in the “Orion Group,” and it 
was considered to be typical of luminosity class iB2. The apparent 
distance moduli of the comparison stars range from 5.5 to 13.3. As 
the mean of the two determinations of the modulus of P Cygni is 
g.3, the selection of comparison stars can be considered satisfactory 
for investigating the amount of selective absorption. 


IV. THE COLOR INDICES 


For the determination of the color indices of the comparison stars 
and of P Cygni, the method of effective wave lengths was used. A 
coarse wire grating having approximately equal bars and spaces 
was used in connection with the 24-inch reflector, which was dia- 
phragmed to 15 inches to improve the photometric field. Eastman 
Ortho Press plates were used. The exposures were made at the same 
altitude as the pole, and the transformation to color indices was 
made by means of direct comparisons with the international color 
indices of the North Polar Sequence. The exposure times were, in 
general, the same on the field of P Cygni and on the pole; they 
ranged in duration from a few seconds to twenty minutes. The ques- 
tion of possible systematic errors depending on the variation of ex- 
posure time was investigated, and it was found that any such errors 
were so small that they were within the mean error of one observa- 
tion. 

The combination of aluminized reflector and orthochromatic emul- 


8 Thid., 47, 463, 1918. 
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sion was found to be very sensitive to slight changes in color index, 
so that the internal accuracy of the color indices is good. A tabu- 
lation of effective wave lengths and measured distances between 
the maxima of the first-order spectra is given in Table 1, where the 

















TABLE 1 
—— as - 
Star | Sp. Cy 2D Neff 
NPS 2 “ome Ao —0.08 402 0.433 
3 Fo +o.18 414 446 
38 F2 +0. 31 420 452 
ir | Ma +1.59 468 0.504 





columns give: (1) the star’s designation in the North Polar Se- 
quence; (2) the Henry Draper spectral type; (3) the color index on 
the International System; (4) the measured distances between the 
maxima of the first-order spectra in w; and (5) the effective wave 
length, also expressed in yw. There is a change in distance of 66 u 
between Ao and Ma; the uncertainty of measurement is about + 2 u. 
Some of the comparison stars were observed only once, while the 
color index of others depends on from two to five observations. 


V. THE INTRINSIC COLORS 

As the North Polar Sequence does not contain stars of early B 
spectral types, the system of normal color indices for such stars is 
rather uncertain. The procedure adopted was to adapt the intrinsic 
colors derived in an earlier paper.’ These colors were on the system 
of Bottlinger’s photoelectric measures. They were first transformed 
to King’s system (King’s photographic magnitude minus the visual 
magnitude in the Revised Harvard Photometry) from a plot of stars 
common to both observers. A scale and zero-point correction was 
then applied to reduce the system as nearly as possible to the North 
Polar Sequence. This was done by making the interval Ao gKo, 
1.12 mag. and assigning bright Ao stars a color index of zero. The 
system of intrinsic color indices so derived is given in Table 2. As 
the bright stars of the Polar Sequence are reddened by about a tenth 
of a magnitude, this correction was applied to the colors derived for 
the comparison stars to determine the amount of interstellar red- 
dening. 


9 Ibid., 87, 471, 1938. 
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VI. THE OBSERVED COLOR EXCESSES 
The resulting color indices and color excesses for the fifteen com- 


parison stars and for P Cygni are given in Table 3. The columns 


TABLE 2 


1} 
| 
| 


Type ¢ Type Cc 





Bo —0.58 dF 5 +0. 37 
B2 | —0.45 dF8 0.48 
Bs -0.30 dGo 0.55 
BS | —= 22 dGs5 } 0.69 
\o ° dKo 0.90 
2 | +0.13 dK5 eee 
Ax... | +o0.21 dK6 1.47 
Fo | +0.23 
gk’5 ca| Spor cF5 0.2 
gk8 | -+0.58 cF8 0.75 
gGo +0.74 cGo bt 
gGs5 +0.92 cGs 1.30 
gKo +1.12 cKo 1.58 
gK2 | +1.26 cK2 iy 
gKk5 +1.67 || cK3 1.99 
gMo... +1.88 cMo.. +2.15 
TABLE 3 
— _ — 
HD Mvyis Sp. M m—M Cy AC 
| 
193184. . 8.5 dF 5 +-3.0 5.5 0.23 —0.04 
193012 Q.2 A2 +1.5 7 a 15 + .12 
228944 9.8 A3 TI.5 6.3 18 + .12 
228623... 10.1 A3 +1.5 8.6 26 + .20 
228793 10.1 dA2 +I. 5 8.6 20 + .17 
193537 9.2 Ao +0.5 8 18 + .28 
228618 10.7 As +2.0 8.7 24 + .13 
228508 9.8 Ao +0.5 9.3 22 Bees v. 
228906 10.1 Ao +0.5 9.6 17 + .27 
228932 10.1 B8 —O.5 10.6 15 + .48 
228807... 8.3 iB2 —3.1 11.4 16 + .71 
193516 8.8 iB2 -3.1 11.9 21 + .76 
193443 7-7 i+Bt 4.9 12.4 i | .77 
193182 a1 gBr —5.5 12.6 18 + .80 
193444 9.0 i+B2 a 13.3 22 + .77 
} 
P Cygni 4.9 Bip | Ho, 0.22 +0.86 























give: (1) the Henry Draper number; (2) the visual magnitude; (3) 
the spectral type (if there is no prefix the type is that of the Henry 
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Draper Catalogue; the other stars were reobserved at Yerkes and were 
assigned to luminosity classes; the prefix ‘i’? corresponds to stars 
of intermediate luminosity, similar to y Orionis; the classification 
“i+” is assigned to stars intermediate in luminosity between y 
Orionis and the supergiant x? Orionis); (4) the absolute magnitude; 
(5) the apparent distance modulus; (6) the color index on the Inter- 
national System; and (7) the color excess, found by comparing (6) 





qT q T ' ' tT ' q q t ' | 


+0.yok 


Color excess 














o~r — ee 
1 i 1 l l iT l l 1 1 at 1 
a $s © 37 © oo 20 3% 12 34 -34. 25 


Apparent distance modulus 


FIG. 1 


with Table 2 and adding o.1 mag. to correct for the selective absorp- 
tion in the Polar Sequence. 

The apparent distance moduli and color excesses are plotted in 
Figure 1. The star P Cygni seems considerably redder than it 
should be if the reddening is due to interstellar selective absorption 
exclusively. If Figure 1 is taken at its face value, we find that inter- 
stellar absorption accounts for a reddening of about three-tenths of a 
magnitude, while’ the residual color excess of about half a magnitude 
must originate in the vicinity of the star itself. 

VII. CONCLUSIONS 


There are, however, several factors which must be considered be- 
fore we accept the quantitative data of Figure 1. They are: (a) the 
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ul 


spectral type of P Cygni; (b) the system on which the color excesses 
are expressed; (c) the possibilities of systematic errors in the spectro- 
scopic luminosities on which the apparent distance moduli are based. 

a) The spectral type of P Cygni.—This is a factor which may be of 
considerable importance. The spectral type is based on the ordi- 
nary strong absorption lines; these lines, however, do not originate 
in the star itself but in an expanding shell where conditions of excita- 
tion may be very different from those in a normal stellar reversing 
layer. If the continuous radiation originates in the star itself 
which seems to be the case—our intrinsic color as determined from 
the spectral type may be considerably in error. We may make the 
problem a little more definite by stating that if the continuous radia- 
tion near the star has a spectral distribution similar to a star of 
spectral type Ao, then there is no reddening which could not be ex- 
plained by interstellar selective absorption. That this could be the 
case is, of course, quite unlikely, and it seems highly probable that 
at least part of the color excess is due to material near the star. 

b) The system on which the color excesses are ex pressed.—A different 
kind of difficulty, but one which complicates the comparison of vari- 
ous determinations of color excesses, is a marked dependence of the 
amount of the color excess on the spectral range observed. The 
values by various observers for five stars are given in Table 4. The 
color indices as observed are listed, and the colors are referred to 
those of stars of the main sequence by means of the comparison stars 
used by each observer. The last horizontal line gives an indication 
of the “‘base line”’ used by showing the interval Ao-gKo. The largest 
interval is that of the spectrophotometric measures, and we find 
that the resulting color excesses are considerably larger than those 
determined by Stebbins and Huffer with the photoelectric cell. On 
the other hand, there seems to be a systematic effect among the 
photoelectric observers themselves which is also dependent on the 
base line used. The stars 55 Cygni and x? Orionis, which are common 
to all observers, show a general decrease in color excess as the base 
line is decreased. The star o Cygni is especially striking; from the 
measures of Becker and Bottlinger it would be considered to be un- 
affected by selective absorption, while the spectrophotometric re- 
sults suggest the opposite conclusion. It is possible that this disa- 
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greement may be due to the same factors which cause the disagree- 
ment between the photographic and photoelectric data of the North 
Polar Sequence. 

c) The spectroscopic luminosities. 
Table 3 are also the most luminous, any systematic error in the abso- 
lute magnitudes of these objects will alter the slope of the distance- 
reddening diagram. It is unfortunate that fainter stars of lower 
luminosity could not be added to the comparison stars to check the 


As the most distant stars in 


TABLE 4 


COLORS OF REDDENED STARS REFERRED TO MAIN 
SEQUENCE INTRINSIC COLORS 
































SPECTRO- STEBBINS (p.-e.) 30TTLINGER 
; BECKER (p.-e.) | 
PHOTOMETRIC* | HUFFER | p.-e.) 
STAR ia a | | 
a ‘or ees 
7 Csp ( | Csp ( | Csp | ( | Csp 
. ames | | | » | | 
P Cygni Bip | 5700°| g3 |+0.06 | ag |—0.08 | a2 | 
55 Cygni cB2 5500 g5 |I+. o5 a8 |— .06 | a3 i=0.54°| a2 
x? Orionis cB2 6100 fo |-0 02: |' jag: j= .07 | a3 58 | ao 
o Cygni cAo | 7700 fo | | — 15 | bo 60 | bo 
v Cephei cA2 | 5700 go | | +0.01 | fo |—0.47 | f2 
— eer 
, | 
Ao-gKo (rex) | 0.74 0.50 | 0.53 
| 





* Miss Payne’s reduction (The Stars of High Luminosity, p. 102) of the results of He asad: Gerasi 
movié, and the Greenwich observers. 
luminosities of the latter. New spectrographic equipment is now 
available which will allow the extension to fainter stars, and we hope 
to strengthen the evidence of the diagram for the more distant stars 
next season. It seems safe to say at the present time, however, that 
a change in slope sufficient to make P Cygni lie on the curve seems 
exceedingly unlikely. 

With due regard to the complications mentioned, we may con- 
clude that the color excess of P Cygni is considerably greater than 
that of other stars in its neighborhood at the same distance from the 
sun, and that a considerable part of the observed reddening is 
probably due to physical conditions in the immediate vicinity of 
the star. 

YERKES OBSERVATORY 
WILLIAMS Bay, WISCONSIN 
January 3, 1939 



































OBSERVATIONS MADE WITH THE NEBULAR 
SPECTROGRAPH OF THE McDONALD 
OBSERVATORY. II* 


OTTO STRUVE AND C. T. ELVEY 


ABSTRACT 

A series of spectrograms in Monoceros and Canis Major reveal a number of new 
emission nebulosities, some of which seem to be unrelated to any individual stars while 
others form distinct groups which are related to O or B stars. A spectrophotometric 
study of Ha and //8 in the outer loop of Orion near BD+1°1126 gives ny = 4 cm-? 
and n; = 3 cm~?. It is estimated that the faintest recorded emission nebulosities have 
three times fewer atoms in the fourth and the third energy-levels of H, per square centi- 
meter. This agrees satisfactorily with former estimates. 

The present series of observations is a continuation of those pub- 
lished in former issues of this Journal.’ The instrument and the 
technique were the same as those used before, but the Agfa Superpan 
Press film was hypersensitized with ammonia, in accordance with 
the method of Bowen and Wyse.’ A marked increase in speed was 
obtained. The numbering of the regions follows chronologically 
that of our previous series. 

The emission features of the spectrum of the night sky were 
relatively faint throughout the present series of observations. This 
was especially valuable in distinguishing the nebular emission line 
(O 11] 3727, which was of even greater use in this series than Ha. We 
are not prepared to state that there is a real systematic difference in 
the ratio |O m|/Ha between the summer Milky Way (Cygnus, 
Cepheus) and the winter Milky Way (Monoceros, Canis Major). 
But there can be no doubt that the large diffuse nebula around 
S Monoceratis is especially strong in |O 1m]. This, incidentally, ex- 
plains our former observation that this nebula is relatively weak 
when photographed in red light.° 

The preponderance of emission spectra in Table 1 is partly due 
to the fact that we have included in our program a number of known 


* Contributions from the McDonald Observatory, University of Texas, No. 9. 


t Ap. J., 89, 119, 1939. 
2 Pub. A.S.P., §0, 305, 1938. 3 Struve, 4p. J., 86, 98, n. 13, 1937. 


517 








518 OTTO STRUVE AND C. T. ELVEY 


diffuse nebulosities for which hitherto no spectra were available. 
The great majority of these extended films of nebulosity in and near 
the Milky Way have emission spectra. However, several regions 
show no trace of nebulosity on the direct photographs, although the 
spectrograph reveals conclusive evidence of emission by interstellar 
gaseous matter. Only one of the regions listed in the table was 
located far from the Milky Way. It showed no emission. In addi- 
tion, the evidence of the comparison spectra is instructive: they 
were usually taken 30°-40° north of the Milky Way region. This 
would place them about 10° north of the eastern edge of the Mi'ky 
Way. Not one of these exposures reveals a trace of emission. 

The principal purpose of the observations was to obtain a de- 
termination of the number of hydrogen atoms per square centimeter 
of a standard nebula. The value used in our former work (1, = 10 
atoms/cm?)4 depended upon a comparison with the Orion nebula for 
which Ambarzumian* had determined the numbers of hydrogen 
atoms in the fourth, fifth, and sixth energy-levels. An independent 
determination is greatly needed because of the astrophysical sig- 
nificance of these quantities.° 

We have chosen for this purpose the outer loop of Orion, near the 
star BD+1°1126. The nebula is relatively strong at this point and 
is fairly uniform over the length of our slit. We estimate that the 
faintest emission nebulae recorded by us are about three times fainter 
than the region selected as standard. The spectral type of the stand- 
ard star is G5, and the photographic magnitude is 6.95, according to 
the Henry Draper Catalogue. Dr. C. K. Seyfert has very kindly de- 
termined the magnitude with the 6-inch UV refractor. His results 
are 

BD + 1°1126 photographic magnitude 6.67 
photovisual magnitude 5.92 
Dr. W. W. Morgan has redetermined the spectral type at the Yerkes 
Observatory and has found it to be G5 on his system. The lumi- 


4 Ibid., 88, 367, 1938. 
5Zs.f. Ap., 6, 112, 1933. 
6 See Struve, Proc. Nat. Acad. Sci., 25, 36, 1939; Struve, Wurm, and Henyey, /bid., 


p. 67; Struve, Zs. f. Ap., 17, 316, 1939. 
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nosity classification by Morgan makes it a subgiant or a dwarf, and 
the color index on the international system should therefore be 
+o.7 mag—in close agreement with the determination by Seyfert. 

The procedure was quite similar to that adopted by Ambarzumi- 
an. The star served as guiding star for the nebula. Hence, there are 
no corrections for extinction. The exposure was obtained on January 
18/19 (see Region No. 41), and the star was trailed uniformly over 
a distance of 11’ along the slit. Sensitometer exposures were used to 
construct the characteristic curves separately for each wave length. 
The width of the slit corresponded to 2‘ 4. 

The only important difference between our work and that of 
Ambarzumian consisted in the fact that in our instrument the slit 
acts as a diaphragm which reduces the intensity of the star but 
leaves the nebula unaffected. For the equivalent slit width of 50 mm 
the correction resulting from this effect corresponds to somewhat 
less than a factor of 2. We shall adopt a factor of 2 in order to com- 
pensate for small losses of starlight at the ends of the trailed image 
and because of optical imperfections in the plane mirrors. 

Although the spectral type is not exactly that of the sun, the re- 
sulting error is insignificant if we use the photographic magnitudes 
of star and sun for HB, etc., and the photovisual magnitudes for 
Ha. 

We apply Ambarzumian’s expression, corrected for the dia- 


phragming effect of the slit: 


En = log P, — 0.4(my — mo) + log S — log 2, 


' 
log Eo, 
where £,, is the energy radiated by 1 cm? of the nebula per second 
in unit solid angle, in the monochromatic line A, Ee, the energy 
radiated by 1 cm? of the sun’s surface per second per angstrom unit 
in unit solid angle, near wave length \, P, the ratio of intensity in 
the monochromatic nebular image, per square minute of arc, to 
intensity in the star spectrum at neighboring wave lengths, per 
angstrom unit, over the full width of the spectrum, expressed in 
minutes of arc, and S the surface of the sun in square minutes of arc. 
Since the linear dispersion on our plates is approximately 1000 
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A/mm at Ha and 440 A/mm at H®, we can easily derive P, for 
both lines and then compute log E,,/ Eo ,. The results are: 


En 

HB: log . a 10.26, 
Eo, 
E, 

Ha: log —* = — 9.56. 
Eon 


For 8 in the Orion nebula, Ambarzumian found log E,,/ Heo, = 
— 7.95. Hence, our nebula is about 2000 times fainter. Remember- 
ing that we have estimated our faintest nebulae to be about three 
times fainter than the standard in the outer loop of Orion, we 
conclude that our former ratio of 1000 was too small by a factor of 6. 
Hence, the number of atoms in the fourth energy-level of 7 for the 


outer loop of Orion is 


my =4cm”’, 


and for the faintest recorded nebulae 


n,= Icm”?. 
For Ha, which was not observed by Ambarzumian, we first compute 
the value of £,, and then derive the number of atoms in the third 
energy-level: 
4mw En 
nN, == Petrik 
hvA,, 


where A,, = 0.43 X 10*°sec'. The result is, for the outer loop of 


3 


Orion, 


This result agrees reasonably with the one used in our previous 
computations (”; = 10 cm”), but the precision of the measurement 
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wal 
nN 
ar 


is, of course, very low. It is, therefore, not disturbing that n, is 
slightly smaller than m,. Both results may easily be in error by a 
factor of as much as 10. 

No. 36.—Guiding star 15 S Mon. 1939, Jan. 14/15. Exp. 3545™. 
Comparison mirror Aé = +30°. Slit eq. 40 mm. [O 11] 3727 is very 
conspicuous; /a is strong; HB, Hy, H6, and N, are present. The 
emission is slightly stronger west of the star than at the eastern end 
of the slit. This nebula has already been observed by Greenstein 
and Henyey.’? The most conspicuous feature is the great strength 
of [O | 3727 and the relative weakness of Ha. The spectral type of 
the guiding star is O7s. 

No. 37.—-Guiding star BD+o° 1376. 1939, Jan. 14/15. Exp. 4°. 
Comparison mirror Aé = +30. Slit eq. 40 mm. This region is 
about 1° east of No. 36, and the spectrum is closely similar. [O 1] 
3727 is very strong; Ha is relatively weak; V, is very weak; 16 and 
Hy are faint. The Balmer decrement seems to be low.’ The emission 
weakens toward the east but is still fairly conspicuous at the eastern 
end of the slit (45’ west of the guiding star). This agrees with the 
extent of the large nebulosity photographed by Ross. 

No. 38.—Guiding star BD+10°1335. 1939, Jan. 15/16. Exp. 
4"4™. Comparison mirror Aé = +30. Slit eq. 40 mm. This region 
is about 4° east of No. 36. There is no definite emission, but we sus- 
pect an exceedingly faint trace of [O 1] 3727. This may, however, be 
the night-sky line of the same wave length. 

No. 39.Guiding star BD+10°1335. 1939, Jan. 15/16. Exp. 
4¢7™. Comparison mirror Aé = +30°. Slit eq. 40 mm. This region 
is about 1 west of No. 36. The small bright nebulosity, NGC 2245, 
which surrounds the guiding star shows no trace of emission, and 
its spectrum must be purely continuous. This agrees with the ob- 
servation by Hubble.’ There is, however, superposed a weak emis- 
sion spectrum, strong at the eastern end of the slit and weak at the 
western end, which evidently represents an extension of the large 
nebula around S Mon. [O u] 3727 is fairly strong, while Ha is very 
weak. The spectral type of the guiding star is B8, according to the 
Henry Draper Catalogue. 

7Ap. J., 87, 80, 1938. 

5 Ibid., 67, Pl. VI, 1928. 9 Tbid., 44, 196, 1916; 56, 173, 405, 1922. 
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No | Object 

36 \Nebula 

| NGC 2264 
37 [Nebula 

| 
38 Milky Way 
30 NGC 2245 
40 IC 448 
41 Outer loop of 

Orion 

42 Milky Way 
43 Milky Way 
44 Milky Way 

| 
45 |Nebula 

fe 2777 

40 
47 
48 





Guiding 


Star 


| 15 S Mon __{1939, 


i1BD+ o 


| 
} 
| 
| 
| 


|BD+10 


|BD+10 
| 


1BD+ 7 


| 


BD+ 1 


| 
| 
| BD+ 8 


BD+15 


| 
| 
| 
| 
| 
| 


BD—1I0 


| p Leonis 


376 


wr 
a 
Jt 


I1S9Q 


1120 


30 7 CMa 


1502] 


| 
| 


Milky Way |BD—10 1892| 


| 


Milky Way |BD— 8°1862| 


Date 


Jan. 


Jan. 


Jan. 


Jan. 


Jan. 


Jan. 


Jan. 


Jan. 


Jan. 


Jan. 


Jan. 


Jan. 
Jan. 


NEBULAR SPECTROGRAPH 


Identification * Result 


14/15|Ross 35; 168 mm_ {Strong emis 
from right; 162 sion 
mm from bottom 


| | 
| 


Strong emis- 


14/15] Ross 35; 148 mm 
sion 


from left; 160 
mm from bottom} 


| 


| 
5/16] Ross 35; 107 mm_ {Probably no 
5; 1 UDI) 
from left; 165 | emission 


mm from bottom| 
| 


5/16] Ross 35; 126 mm__|Emission 
from right; 162 | 
mm from bottom} 


16/17) Ross 35; 139 mm Emission 
from right; 122 
mm from bottom} 
18} Ross 34; 114 mm__|Emission 
| from left; 128 | 
| mm from top 


| 
| 
| 
} 


4 


| | 
17/18|Ross 35; 167 mm__|Strong emis- 


from right; 135 | sion 
| mm from bottom} 
i : 
19/20] Ross 35344 mm__|Strong emis- 
from right; 79 sion 
| mm from top 
19/20] Ross 38: 112 mm__|Emission 
| from left; 146 
| mm from bottom) 


| 
| 
| 


20/21|/Ross 37; 82 mm Strong emis 
| fromright; 164 | sion 

| mm from top | 

} 

| 

| 


22/23| Ross 37; 95 mm Emission 
| from right; 159 _ | 
| mmfromtop — | 
| 


| 
29/23)... ; No emission 
| | 
| a) 
23/24| Ross 37; 139 mm__|Emission 
from right; 134 
mm from top 





* The co-ordinates are only approximate. Differences of the order of 1 mm in each co-ordinate are pos- 


sible on different prints of the same photograph of Ross’s Aflas of the Milky Way. 
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No. go.—Guiding star BD+7°1337. 1939, Jan. 16/17. Exp. 5°. 
Comparison mirror Aé = +30”. Slit eq. 40 mm. This region is 
about 3° southwest of No. 36. [O 11] 3727 is a very strong emission 
line, but Ha is rather weak. The oxygen line is of uniform bright- 
ness along the entire slit and must belong to the large nebula around 
S Mon. Since there is no strengthening of the emission lines near 
the guiding star, we believe that the faint nebula, IC 448, sur- 
rounding it has a continuous spectrum. The spectrum of the guiding 
star is Aop and the lines are described in the Henry Draper Catalogue 
as very narrow and sharply defined, resembling in intensity those of 
n Leonis. 

No. 41.—Guiding star BD+1°1126. 1939, Jan. 17/18. Exp. 
4¢5™. Comparison mirror Aé = +40. Slit eq. 40 mm. 1939, Jan. 
18/19. Exp. 45. Slit eq. 50 mm. Ha is very strong; [O m1] 3727 is 
somewhat fainter, confirming the reality of the unusual strength of 
the oxygen line (or weakness of Ha) in the large nebula around 
S Mon. HB and Hy are faint; N, is not visible on our plate. This 
observation confirms in every detail that by Greenstein and Hen- 
yey’’ which was made slightly farther south. 

No. 42.—Guiding star BD+8°1438. 1939, Jan. 17/18. Exp. 
45™. Comparison mirror Aé = +30°. Slit eq. 40 mm. [On] 3727 
is exceptionally strong; Ha is fairly strong; H® and Hy are weak. 
The lines are fairly uniform along the slit. This is an extension of the 
nebula around S Mon. 

No. 43.—Guiding star BD+15°1075. 1939, Jan. 19/20. Exp. 
3"30™. Comparison mirror Ad = +30°. Slit eq. 40 mm. This re- 
gion was chosen because of the strong emission nebulosity revealed 
in red light, by H. A. Lower." This fairly large, round nebula seems 
not to have been recognized previously, but a careful inspection of 
Ross’s plates No. 35 and 33 shows an exceedingly faint trace of it. 
Barnard’s photograph” does not show it. The nebula must be very 
faint in blue and violet light. This is substantiated by our spectrum: 
Ha is very strong in emission; [O 11] 3727 is relatively weak; H6 and 
Hy are very faint; there is no trace of N,. The emission lines are 
strongly concentrated toward the middle of the slit, suggesting that 

10 Thid., 87, 79, 1938. 

1 Tbid., 89, Pl. XIV, 1939. 12 Pub. Lick Obs., 11, Pl. 25, 1913. 
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the brightest part of the nebula has a diameter of 15’-20’. However, 
Ha extends to both ends of the slit. The nebula is considerably 
larger than should be expected from Hubble’s relation. The spec- 
trum of the guiding star is not known. The nebula may be asso- 
ciated with the slightly brighter star, BD +15°1079 = HD 41997, the 
spectrum of which is given as B in the Henry Draper Catalogue. 
The absorption lines are described as showing slight contrast to 
other portions of the spectrum. The photographic magnitude is 8.3. 
This interesting nebula deserves further study. 

No. 44.—Guiding star 307 Can. Maj. 1939, Jan. 19/20. Exp. 
3"52™. Comparison mirror Aé = +50°. Slit eq. 40mm. This region 
contains no known nebulosity, but a slight milkiness on Ross’s 
photograph suggests the possible existence of a very large diffuse 
nebulosity which merges into a large obscured region on the east. 
Ha is strong in emission and so is [O n| 3727; HB and Hy are 
faint; V, is not definitely seen but may be present. The lines are 
slightly stronger near the eastern end of the slit than in the im- 
mediate vicinity of the guiding star. The spectral type of 7 Can. 
Maj. is Oo, and that of the neighboring star, 29 Can. Maj., is O7e. 

No. 45.—Guiding star BD—10°1862. 1939, Jan. 20/21. Exp. 
3"54™. Comparison mirror Aé = + 40°. Slit eq. 40 mm. Ha and 
[O 11] 3727 are very strong and decrease slightly from the guiding 
star toward the eastern end of the slit; HB, Hy, and H6 show the 
same structure; N, is clearly visible, but its structure differs from 
that of the other lines; it is absent west of the guiding star, is faint 
immediately east of it, and increases slightly in intensity toward the 
eastern end of the slit, 45’ from the star. The spectrum of the guid- 
ing star is B3. The strengthening of V, toward the east suggests that 
the excitation may come from the O7s star™ BD—10°1892 = 
HD 54662. 

No. 46.—Guiding star BD—10°1892. 1939, Jan. 22/23. Exp. 
2290™, Comparison mirror Ad = + 40°. Slit eq. 40 mm. The Ross 
Atlas shows no trace of nebulosity in this region, although it lies 
only about 1° east of the diffuse nebula IC 2177 (Region 45). Ha 
and [O u| 3727 are strong over the entire length of the slit, weakening 
slightly toward the east; NV, is equal to H@, thus confirming the 


13 C. H. Payne, The Stars of High Luminosity, p. 64, 1930. 
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result obtained for Region 45. The spectrum of the guiding star 
is O7s. 

No. 47.—Guiding star p Leonis. 1939, Jan. 22/23. Exp. 2"40™. 
Comparison mirror Aé = +30°. Slit eq. 40 mm. This Bop star has 
a galactic latitude of +54°. There is no emission. 

No. 48.—Guiding star BD—8°1862. 1939, Jan. 23/24. Exp. 
3h5™. Comparison mirror Aé = + 40°. Slit eq. 40 mm. This region 
is about 4° northeast of No. 46 and is located in a dense region of the 
Milky Way, without any indication of nebulosity on the direct 
photograph by Ross. [O 11] 3727 is present as a faint but well-de- 
fined emission line which decreases slightly in strength near the 
eastern edge. There is only a suspicion of bright Ha, but this plate 
is slightly out of focus in the red part of the spectrum. The presence 
of emission in this normal region of the Milky Way is perhaps even 
more startling than the presence of emission in the star clouds of 
Cygnus and Cepheus."4 

McDonaLpD OBSERVATORY 


AND 
YERKES OBSERVATORY 


February 2, 1939 


"4 Dr. W. H. Wright has called our attention to the fact that the strong Ha line coin- 
cides on spectrograms of small dispersion with the forbidden doublet of NV 1, \ 6548 and 
\ 6584, and there is good reason to suspect that the red images of certain nebulae are due 
principally to this doublet (Pub. Lick Obs., 13, 217, 1918). Since the other Balmer lines, 
HB, Hy, etc., are also often present on our spectrograms, Ha is doubtless a fairly strong 
emission line; but it is entirely possible that a large fraction of the light comes from the 
nitrogen doublet. If this conclusion is correct, the number of atoms which we have 
derived for the third quantum level of hydrogen must be somewhat reduced. However, 
it will not be altered sufficiently to require a revision of the computations because we are 
here concerned only with a precision corresponding to a factor of 3 or 4. Dr. Wright’s 
suggestion probably accounts for the remarkable changes in the relative intensities of 
Ha and of other Balmer lines which we had heretofore interpreted as changes in Balmer 
decrement. Stoy (Lick Obs. Bull., 17, 181, 1935) assigns slightly weaker intensities to 
the [.V 11] lines than to Ha, in the nebulae observed by him, but it is quite possible that 
the combined intensity of the two [.V 1] lines in our tenuous nebulae is of the same 


order as that of Ha. 
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ABSTRACT 


The discovery, by Struve and Elvey, of extended areas in the Milky Way in which 
the Balmer lines are observed in emission suggests that hydrogen exists, in the ionized 
state, in large regions of space. The problem of the ionization and excitation of hydro 
gen is first considered in a general way. An attempt is then made to arrive at a picture 
of the actual physical state of interstellar hydrogen. It is found that the Balmer-line 
emission should be limited to certain rather sharply bounded regions in space surround 
ing O-type stars or clusters of O-type stars. Such regions may have diameters of about 
200 parsecs, which is in general agreement with the observations. Certain aspects of the 
problem of the ionization of other elements and of the problem of the relative abundance 
of the elements in interstellar space are briefly discussed. The interstellar density of 


hydrogen is of the order of NV = 3 cm~3. The extent of the emission regions at right 
angles to the galactic plane is discussed and is found to be small. 


The recent discovery, by Struve and Elvey,' of extended regions 
in the Milky Way showing hydrogen-line emission has opened up 
new and highly important possibilities for the study of the properties 
of interstellar matter. From the observed strength of Ha in the 
emission regions, Struve? has calculated the density of interstellar 
hydrogen. Also, he has analyzed the problem of the ionization of 
interstellar calcium and sodium, taking account of the presence of 
ionized hydrogen. 

In the present paper the problem of the ionization and excitation 
of interstellar hydrogen is first considered in a general way. Based 
upon the observed intensities of interstellar hydrogen emission, an 
attempt is then made to arrive at a picture of the actual physical 
state of interstellar hydrogen. Finally, certain aspects of the prob- 
lem of the ionization of other elements and of the problem of the 
relative abundance of the elements in interstellar space are briefly 
discussed. 

II 

Eddington has considered the problem of the ionization of inter- 
stellar hydrogen and has expressed the opinion that, in a normal re- 

™ Ap. J., 88, 364, 1938; 89, 119, 1939. 2 Proc. Nat. Acad. Sci., 25, 36, 1939. 


3 M.N., 95, 2, 1934; Observatory, 60, 99, 1937. 
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gion of interstellar space, hydrogen is entirely un-ionized, the reason 
being that the ionizing ultraviolet radiation is strongly absorbed by 
the interstellar hydrogen. The result of the analysis given below 
tends to confirm Eddington’s view with the modification, however, 
that high-temperature stars, and especially clusters of such stars, 
are capable of ionizing interstellar hydrogen in regions large enough 
to be of importance in problems of interstellar space. 

In the immediate neighborhood of a star, interstellar hydrogen 
will be ionized. With increasing distance from the star the propor- 
tion of neutral hydrogen atoms increases, and hence the absorption 
of the ionizing radiation increases. Ultimately, the ionizing radia- 
tion is so much reduced that the interstellar hydrogen is un-ionized. 
Our problem is to derive an expression for the extent of the ionized 
region as a function of the temperature and absolute magnitude of 
the star and the density of interstellar hydrogen. A somewhat simi- 
lar problem has been considered by the author in another connec- 
tion.4 

Consider a point in interstellar space at the distance s from a 
star of temperature T and radius R. Let the number of neutral hy- 
drogen atoms per unit volume be NV’; the number of hydrogen ions, 
N”; and that of free electrons, V.. The degree of ionization at s 
is determined? by the equation 


NN, _ (2em.)3? 29” ,,,. fe a 7 
N’ =e ! me) 7 (kT)3/2e? veal Ve ‘me "ea, (1) 








Here J is the ionization potential; g’’ and gq’ are the statistical 
weights of the ion and the ground state of the atom, respectively; 
Tis T is a correction factor to allow for the difference between the 
temperature T of the exciting star and the electron temperature 
T., at s (cf. Rosseland’); while e~’“ measures the reduction in the 
ionizing ultraviolet radiation due to absorption, 7, being the optical 
depth from the ionizing star to the point s. Finally, w is the dilution 
factor at s, given by 
w= « : (2) 
4G. P. Kuiper, O. Struve, B. Strémgren, Ap. J., 86, 570, 1937; sec. III, especially 
pp. 593 f. and 612. 


5 Theoretical Astrophysics, chap. xxii, Oxford, 1936. 
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We shall assume that practically all the free electrons are fur- 
nished by ionization of hydrogen, so that NV” = N,. The validity 
of this assumption for actual interstellar space will be discussed in 
section V. Further, let 


N =N’+N", 


N” =<xN, | 

(3 
N’ = (i — x)N > | 
N,. =xN, 


so that x is the degree of ionization of hydrogen. Introducing nu- 
merical values, equation (1) can now be written in the form 


x . I 
N= See (4) 
with 
2 47 T. : 
C, = 19°0-"F =f. we r ee, 
q 1 
° (5) 
5040 
§= 2 
7 


The numerical factor of C, corresponds to the following choice of 
units: 1 parsec = 3.08 + 10"* cm for s; the solar radius for R; and 
cm~$ for NV. 

Let the absorption coefficient for the ionizing radiation per neu- 
tral hydrogen atom be a,. We simplify the problem by assuming 
a, to be independent of the frequency and equal to its value at the 
absorption edge. The relevant range of frequency is relatively 
small, so that the accuracy will be sufficient for our present purpose. 
Then, by definition, 


dt, = (1 — x)Nay + 3.08 © 10% ds , (6) 
the factor 3.08 - 10" being derived from the choice of the parsec as 


the unit of s. 
Equations (4) and (6) define the solution of the problem of deriv- 
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ing the degree of ionization as a function of the distance from the 
star. Solving (6) for 1 — x and substituting in (4), we get 


e "udr, = ron x75? + 3.08 - 108a,uds . (7) 
I 
We introduce the following new variables: 
y= eu (2.339 >a), @& 
N? s ‘ : 
dz = i 3.08 * 10a, ° s7ds, (9) 
I 


with s = o for s = 0, so that 


3C; ) : 1/3 
s= | ——_——_—— o gis, IO 
(x ° 3.08 + 10a, ag 


Then (4) and (7) can be written as 


dy ; 
a oe rH 
dz ; (11) 

I—-x I 
ee SS Og 3? ee EZ 
x? y 12) 
with 

re) aie , 
elk fy een oe ot a 13 
(x 1° (3.08 + Th (13) 


Since 7, = o for s = o, we have, according to (8) and (10), 
y=1 for z=o0. (14) 


In the cases of actual interest, a is a small quantity. When 
a <1,it follows from (12) that, as long as y is not small compared 
to 1,1 — « <1, so that x is equal to 1, very nearly. In that case 
it follows from (11) and (14) that 


y=1-2 for 1—-xXK1. . (15) 
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Consequently, when a is small, ionization is strong until z is nearly 
equal to 1, so that (cf. [15]) y becomes a small quantity. With in- 
creasing z the degree of ionization now decreases very abruptly, so 
that, for z slightly greater than 1, hydrogen is practically un-ionized. 

The value of s corresponding to z = 1, which we shall call s,, is, 


according to (10), 





5, = ( 3Cx ae (16 
sie N? + 3.08 + 10a, ue) 


The result of the analysis can be stated as follows:, for small a the 
ionization is nearly complete up to the distance from the ionizing 














TABLE 1 
S/So 
x? 
| 
a=o.oo!I a=0.0I a=o.!I 
E20:.. 0.000 0.000 0.00 
6.8: I .000 0.988 0.82 
620... I.002 I .009 0.97 
O.#:.. 1.003 1.020 1.05 
0.2 I .004 1.028 ee &: 








star s,, given by (16), while there is almost no ionization for dis- 
tances greater than 55. 

The abrupt change of the ionization can be interpreted in the fol- 
lowing way. Once the proportion of neutral atoms begins to increase, 
the absorption of the ionizing radiation increases, leading to an ac- 
celerated increase of neutral atoms. (For a somewhat more detailed 
discussion of this phenomenon see the investigation by the author 
quoted above.‘) 

The exact dependence of the degree of ionization x upon the dis- 
tance s from the ionizing star has been derived for three different 
values of a by numerical integration according to (11), (12), and 
(14). The results are shown in Table 1. For a = 0.01 the results 
are given in greater detail in Table 2. Finally, it may be noticed 
that for small a the following relation, obtained by integrating (11), 
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(12), and (14), putting = 1, holds very nearly in the transition re- 
gion between almost complete ionization and negligible ionization: 


( x , 
z = Const — «{ ——— + aln — . (a<1,|1—2|K1). (17) 


Numerical values of s + constant, normalized to o for a? = 0.5 
and calculated, according to (17), as a function of x?, are shown in 
Table 3. The table shows, for instance, that with a = o.oo1 the 





























TABLE 2 
a=0.0I 
F4 5/So y=e""u Ty 2 x 
0.0 0.00 I.00 ©.00 ©.000 I.00 
0.2 0.58 0.80 0.22 004 °.99 
0.4 0.74 0.60 0.50 009 0.98 
0.6 0.84 0.41 0.89 O17 0.97 
0.8 0.93 0.22 E82 .037 0.93 
0.9 0.97 0.13 2.06 064 °.88 
1.0 1.00 0.040 3.07 15 O.72 
Bakes 1.03 ©.0018 6.3 0.67 O.II 
TABLE 3 
z+Constant z+Constant 
x? Normalized to x? Normalized to 
o for x?=0.5 o for x2=0.5 
Gebocn se isses. ae od - ==b 76 
22 3.8 7 4.2 
rae 2.6 8 8.6 
oan +1.3 a) i = 26n3 
©.5:.. omre) 


change from x? = 0.5 to x? = o.1 takes place for a change in z equal 
to 0.0053, corresponding to a change in s of about one-sixth of 1 per 
cent. 

For small a the relative extent of the transition region is propor- 
tional to a. From (13) and (16) it follows that as, is independent of 
C,, 1.e., of the properties of the ionizing star. Therefore, as long as 
a is small, the absolute width of the transition region does not vary 
with the properties of the ionizing star. 








1 
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Combining equations (5) and (16), we find 


2q”’ om 
log so = —6.17 + 3 log a) —iloga, — } 6l 
7 3 a 7 3 10g 3 (18) 


+ flog 7+ ZlogR — 2 log N. 


This relation holds for any element. Introducing the proper numeri- 
cal values for hydrogen, viz., g’’/q’ = 3, I = 13.53 volts, and a, = 
6.3 - 10-7 cm-?, we get 


log s5, = —0.44 + 3 log (./72) — 4.510 + 4 tn % 
(19 
+ 2logR — ZlogN .| 


Table 4 gives log s, for hydrogen for R = 1 and N = 1 asa func- 
tion of 7. The temperatures have been so chosen as to correspond 
to the spectral types from Os to Bs, according to the temperature 











TABLE 4 
nee 
| | log So log a 
Sp | log 7 6 for R=1 So | for R=1 

| | and V=1 |} and N=1 
Os | 4.90 0.063 eye 54 parsecs X R?/3N-? 7.40—10 
O06 | 4.80 | 079 1.60 40 | 7.50 
O7 | 4.70 | .100 1.46 29 . 298 
O8 | 4.60 |. .126 I. 20 20 } 7-90 
Oo | 4.50 | .158 1.10 13 | 8.09 
Bo | 4.40 | .200 0.86 7.2 $.22 
Bi... | 4.36 | 210 0.75 5.6 | 8.44 
ae Aste 245 0.62 Aue ) 3Bee9 
B3 A | 4.27 | . 269 ©.49 Z.7 | 8.70 
a oe | 4.23 | .205 0.35 +2 | 8.84 
Sea | 4.19 | 0.324 0.20 | 1.6 | 8.99 

| | | | 





scale recently derived by Kuiper.® In making the calculation the 
factor T.,/T has been put equal to 1. With 74/T = } (cf. Rosse- 
lands), all values of log s, would have to be decreased by 0.10, which 
would correspond to a decrease of the s,-values of 21 per cent. 
Table 4 also gives a according to (13) for R = 1 and N = 13a is pro- 
portional to R-?’3 and to N-"3, 

For any given temperature of the ionizing star, s, increases with 


6 Ap. J., 88, 429, 1938. 
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the stellar radius Ras R?’3. If a region of interstellar space is ionized 
by a cluster of m similar stars, close together, then s, has to be cal- 
culated with an equivalent R equal to n'’?R. This follows immediately 
from equation (2) for the dilution factor w. Consequently, s, for such 
a cluster is equal to 2'/3 times s,, calculated for the individual star in 
the cluster. This may also be expressed by saying that the volume 
of interstellar space ionized by a cluster of stars close together is 
equal to the sum of the individual volumes that would be ionized by 
the stars if placed so far apart in interstellar space that the volumes 
did not overlap. 

Table 5 gives s, for main-sequence stars of spectral types O5-Bs. 
The assumed visual magnitudes’ from which R was computed, using 


TABLE 5 
Sp T | Mets So 

Oe. . 79,0007 | —4™2 | 140 parsecsX N-2/3 
OOo... ; 63,000 —4.1 | 110 

O7 4 50,000 —4.0 87 
ee soa 40,000 =f. 66 

i) 5 32,000 3.60 40 
ee 25,000 — 3.1 26 

Bi Sete oer ds evn 23,000 —2.5 17 
B2... voy 20,000 —1.8 Il 

B3 x eee 18,600 —1.2 | ae 

BA a icBdcwacuseae Teewe —1I.0 | ee 
ae us ash 28.600 —o.8 | 3.7 

IRE cco dea ccs es 10,700 +o0.9 0.5 

















Kuiper’s bolometric corrections,° are also given. The values of R 
range from about 4 to 7 solar radii. An increase in brightness of 
1 mag. would lead to an increase of s, by a factor of 1.36. 

The increase in the ionized volume as one passes from low-tem- 
perature stars to high-temperature stars is so pronounced that it 
may be concluded that the total volume ionized by all high-tempera- 
ture stars (hotter than about B2) is much larger than that ionized 
by all low-temperature stars, in spite of the much greater number of 
the latter. For instance, one O7 star ionizes a volume equal to that 
ionized by about two thousand B3 stars or by about five million Ao 
Stars. 


7R. Trumpler, Lick Obs. Bull., No. 420, 1930. 
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In deriving the values given above, certain simplifying assump- 
VT.4/T has already been discussed. A resulting decrease of the 
So-values on this account by more than, say, 20-30 per cent is not 
very probable. It has also been pointed out that a certain mean 
value of a,, taken over the relevant frequency range of the ionizing 
radiation, should have been used instead of a, at the absorption 
edge. It is estimated that this effect will lead to an increase of s, of 
probably less than, say, 20-30 per cent. Finally, the ionizing radia- 
tion emitted by the interstellar hydrogen has been neglected in cal- 
culating the ionization. An upper limit to the change of s, resulting 
from this effect is obtained in the following manner. Upon absorp- 
tion of a quantum of ionizing radiation by an interstellar hydrogen 
atom there may follow reradiation of a quantum of the same fre- 
quency through a capture process to the ground state. In the case of 
hydrogen under the general conditions in question it is known that 
about one-third of the captures of free electrons by ions are to the 
ground state, the rest being to excited states. Therefore, about two- 
thirds of the absorption processes from the ground state are effective 
in removing entirely the quantum of ionizing radiation. The re- 
maining one-third of the absorption processes also act to reduce the 
resulting ionization, through the scattering involved. The upper 
limit referred to is obtained by neglecting the effect of the scattering, 
thus introducing an effective absorption coefficient equal to } a, in 
(16). The maximum influence is therefore an increase of about 15 
per cent in s,. Altogether, it is estimated that the s,-values calcu- 
lated are sufficiently accurate for a first survey of the problem. It 
would not, however, present very great difficulties to develop a the- 
ory that would take the effects mentioned more accurately into ac- 
count. 

We have tacitly assumed that ionization by electron collisions is 
unimportant compared with photo-ionization. The justification for 
this assumption lies in the fact that, in interstellar space, the radia- 
tion from the ionizing stars is the primary energy source. Without 
making detailed calculations, we assume that any tendency for col- 
lisional ionization to become of importance is checked by a corre- 
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sponding decrease of the electron temperature. The mean free path 
of the electrons is so small as not to interfere appreciably with such 
an adjustment of the electron temperature. 

Applications of the theory developed in this section are given in 
sections IV and V. 


Iil 


We shall consider in this section the problem of the excitation of 
the higher energy-levels of hydrogen in interstellar space. We shall 
limit ourselves mainly to the consideration of the level m = 3, which 
determines the emission of Ha radiation. The discussion is naturally 
divided into two parts, dealing with ionized and nonionized regions, 
respectively. 

In the almost completely ionized region surrounding a high-tem- 
perature star or a cluster of such stars, the mechanisms of excitation 
of one of the higher levels—say n = 3—can be classified as follows, 
leaving out of consideration mechanisms that are entirely negligible: 
(a) capture of a free electron by an ion to the state considered; (6) 
capture of a free electron by an ion to a higher state, followed by 
one or more transitions leading through the state considered; (c) 
Lyman-line absorption from the ground state to the state considered, 
and Lyman-line absorption from the ground state to higher states, 
followed by one or more transitions leading through the state con- 
sidered; (d) Balmer-line absorption from the state m = 2, either di- 
rectly to the state considered or to a higher state followed by cascad- 
ing, as before; and (e) excitation by electron collision from the 
ground state to the state considered or to higher bound states, fol- 
lowed by cascading. 

In investigating these mechanisms we can largely make use of the 
results of investigations of the analogous problem in the case of 
planetary nebulae. 

Consider first mechanism (a). The effective cross-section 7o?(m) 
for capture of a free electron of kinetic energy E, to the state with 
quantum number of hydrogen is given by a well-known formula 
(cf., e.g., Cillié’) which follows immediately from the Sugiura- 


* M.N., 92, 820, 1932; 96, 771, 1936. 
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Gaunt formula for the continuous absorption cross-section and the 
Milne relation between continuous absorption cross-section and cap- 
ture cross-section for reverse processes: 


I 
2(m) = 1. ° —33 & Vo —. 20 
o*(n) a nv E, (ae) 


where g is the Gaunt correction factor, which is approximately equal 
to 1 (g = 0.91 for nm = 3), while », and v are the frequency of the 
absorption edge and of the radiation emitted at the capture, respec- 
tively (so that hv = hy, + xn, when x, is the binding energy of 
the state considered). 

With the aid of the well-known gas-kinetic formula, giving the 
number of captures corresponding to a given capture cross-section 
(cf., e.g., Fowler’), and using the mean lifetimes for the various 
states of hydrogen, as tabulated, for example, by Bethe,’? the num- 
ber of atoms per unit volume J, in any state can now be calcu- 
lated easily (cf. Cillié*). For 2 = 3 one finds 


pe = 6.3° 107+ T-/?+ J (excitation by mechanism [a] only). (21) 
Here J denotes the proper mean value of the factor »,/v~ in (20). For 
very small temperatures, J is approximately equal to 1; for T = 
sooo’, J = 0.81; for T = 20,000°, J = 0.57; and for T = 50,000", 
J = 0.40. 

It may be noted that, in the particularly simple case of excita- 
tion by mechanism (a) there is no difficulty at all in treating the 
substates of the state with m = 3 separately. When this is done, the 
numerical factor in (21) is increased by a factor of about 1.5, the 
number JN, being so defined that the actual number of transitions to 
the state m = 2 is given by the product of JN, and the standard tran- 
sition probability 7;,, obtained by averaging according to the statis- 
tical weights. 

The case of excitation by both mechanisms (a) and (6) has been 

9 Statistical Mechanics, Cambridge, England, 1936. 


10 Handb. d. Phys., 24, 1, chap. lii, 1933. 
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treated by Menzel and by Menzel and Baker." These authors give 
the ratio V,,/N’’N, in terms of the corresponding ratio for thermo- 
dynamical equilibrium. In thermodynamical equilibrium we have 


N, rr . *“¥* . 
yy. = 10 14-44+1.508 6 T—3/2 (thermodynamical equilibrium). (22) 


In the case of excitation by mechanisms (a) and (6) 


N por ‘esti 
Vy, = 54) * Io 14:44 1.508 , T 3/2 } 
sl > (23) 


(excitation by mechanisms [a] and [6}), 


where b{*? is tabulated by Baker and Menzel" (their hypothesis A,) 
as a function of the electron temperature 7.;. The values of , given 
by (23) are naturally larger than those given by (21). For Ty = 
5000° the increase is given by a factor about equal to 2; with increas- 
ing temperature the factor tends toward 1. 

The case of excitation by mechanisms (a), (6), and (c) has been 
considered by Cillié® and by Menzel and Baker." Following Zan- 
stra,’ they treated the case in which the exciting Lyman-line radia- 
tion, as well as the ionizing radiation, which determines the relative 
number of atoms in the ground state (from which the excitation con- 
sidered takes place), is given by a common dilution factor times the 
Planck intensities for a common temperature. In this case the tran- 
sitions from any excited state to the ground state are exactly bal- 
anced by the reverse Lyman-line excitations, and the treatment is 
analogous to the treatment of the case (a) + (b), except that all 
transition probabilities to the ground state are put equal to zero. 

In the particular case of excitation by mechanisms (a), (b), and 
(c) just considered, 


- = O(B) 6 youtastt509 6 [3/2 
(24) 


(particular case of excitation by mechanisms [a],[6], and [c]), | 


1" T). H. Menzel, Ap. J., 85, 330, 1937; D. H. Menzel and J. G. Baker, Ap. J., 86, 70, 
1937; J. G. Baker and D. H. Menzel, Ap. J., 88, 52, 1938. 


12 Pub. Dom. Ap. Obs., Victoria, 4, 200, 1931; Zs. f. Ap., 2, 1, 1931. 
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where b{) is tabulated by Baker and Menzel (their hypothesis B) 
as a function of the electron temperature. 

With regard to the importance of excitation by mechanism (d), 
we shall only consider two limiting cases. If the excitation of the 
state m = 2 is normal for interstellar space, the role of mechanism d 
is entirely negligible. If, however, the excitation of the state n = 3 
is supernormal to such an extent that the ionized region becomes 
opaque in the frequencies of the Balmer lines, then the relative popu- 
lation of the state m = 3, excited by mechanisms (a), (0), (c), and 
(d), approaches the value given by expression (22), valid for thermo- 
dynamical equilibrium, with the temperature about equal to that of 
the ionizing star. (For the possibility of superexcitation of the state 
n = 2, owing to the metastability of the 2s-substate, see Chandra- 
sekhar and Stoy, as quoted by Cillié.*) 

With regard, finally, to the importance of excitation by mecha- 
nism (e), it is easily seen that in an almost completely ionized region, 
and for temperatures of the order of those to be expected in actual 
applications to problems of interstellar space, the mechanism is un- 
important. The average capture cross-section for the capture of a 
free electron by an ion, to the state m = 3, will, for the temperatures 
in question, be of the order of magnitude 10". The cross-section for 
excitation by an electron of a neutral atom in the ground state to 
the state m = 3 is of the order of magnitude 10~”, if the energy of the 
electron is high enough to excite the state m = 3 at all. The fraction 
of electrons that have sufficient energy for excitation is of the order 
of magnitude 10-3. When, as in the case under consideration, the 
ionization is so high that the number of ions is, say, 107-104 times 
the number of neutral atoms, then mechanism (e) will be unimpor- 
tant compared with mechanisms (a) and (0). 

It is interesting to note, as can easily be shown, that even in such 
cases of less complete ionization, when mechanism (e) is of impor- 
tance, the ratio V,/N’’N. will approach the value given by equation 
(23), the reason being that, for high-level states like m = 3, the cross- 
section for electron ionization from the ground state is considerably 
larger than the cross-section for electron excitation of that particu- 
lar state from the ground state (cf. Bethe’). In interstellar space, 
moreover, we assume that electron ionization is unimportant com- 
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pared with photo-ionization as a consequence of an adjustment 
of electron temperature. This effect also limits the importance of 
mechanism (e) in the case considered. 

Table 6 illustrates equations (21), (22), (23), and (24), valid for 
the various cases considered. However, the table does not give the 
ratio V,/N”N, directly. It gives, with a view to the application in 
the following section, the value of V,, calculated with V’”’ = N, and 
N,= 3-10. For other values of V,, NV. may be obtained by 


TABLE 6* 























| tite ? 
| (1) | (2) | 3 (4) (4 ) 
| Thermody 
namical eee 
Eauilit | Special Excitation 
Squill ri . ° ‘ . 
e T : S Case of | Excitation Excitation (a); Sub- 
“- igi | Excitation a)+(6) (a) states Con- 
ase Oo . P . . 
acts % | (a (b)+(c)| (Eq. [23]) (Eq. [21]) sidered 
“xcitation > | 
(Ex 24]) | Separately 
(a) +(6)+ | | I ’ 
| c)+(d) | | 
| | ' 
1.0 ,O40 o.10cm | r.ocm'3; 1.4cm3 2.0cm 3 I1.7cm 3 
| 
0.8 | 6,300 0.16 a3 | ¥.6 / ey - 1.5 
0.6 | 8,400 0.28 | Sit ee 2.4 2:6 
0.4 12,600 0.55 } 1.4 } 2.2 2.5 2.3 
0.2 25,200 ae | 2.0 | 3.0 3.8 3.1 
O.1 | 50,400 2.5 | 2.5 | 4.3 5.2 4.2 
| 
* The table gives VN, for N’’=N, and N3=3 ° 107?! cm73, 


multiplication with the square root of the factor by which J, differs 
from 3 - 10°”. 

It is interesting to note that for the whole range of excitation 
mechanisms considered, and also for the relevant range of tempera- 
tures, the values of V, come out approximately the same. (As was 
noted above, in the special case of excitation [a] + [b] + [c] + [d] 
the temperature approaches the temperature of the ionizing star, so 
that the relatively small numbers in the upper left corner of Table 
6 are irrelevant. For mechanism [a] + [b], say, the temperature in 
question is the electron temperature, and so may be considerably 
lower than that of the star; but here the change in the tempera- 
ture affects V. very little.) The relative insensitivity of the ratio 
N,/N’N, here encountered is not surprising. We may mention the 
analogous example of the approach of the central intensities of stellar 
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absorption lines to the value for thermodynamical equilibrium, for 
temperatures high compared with the binding energy. 

Further, it is of interest to note that for mechanisms (a), (0), (c), 
and their combinations the ratio V,/N’’N, is independent of the dis- 
tance from the ionizing star. This is not true for mechanism (e), 
which, however, as we have seen, is unimportant in the present con- 
nection. 

We shall now consider briefly the actual mechanism of excitation 
in the ionized regions around the high-temperature stars. The prob- 
lem is not quite simple, but it is suggested that, in the major part of 
the ionized region, Lyman-line radiation is present almost only 
in the form of La, the radiation in the other Lyman lines having been 
converted by repeated fluorescence processes into La, Balmer lines, 
Paschen lines, etc.; also, that the ionized region is transparent in the 
frequencies of the Balmer lines (cf. above), so that the excitation of 
the state m = 3 and of higher states is by mechanism (a) + (d) 
(Table 6, case 3). 

In the paper quoted at the beginning of this article, Struve? cal- 
culated the ratio N,/N’’N, in an ionized region of interstellar space 
by a procedure leading to equation (22) (Table 6, case 1), namely, 
the use of the Boitzmann and the Saha formulae with a common 
dilution factor which drops out in forming N,/N’’N,. It follows 
from the previous discussion that the value of the ratio obtained by 
Struve will not differ very much from those obtained on any of the 
other hypotheses discussed. The improved method for calculating 
N,/N”N, briefly outlined by Struve is essentially equivalent to that 
used here in case 3, i.e., the case which probably corresponds most 
closely to actual conditions. 

We finally have to discuss the excitation of hydrogen in the prac- 
tically un-ionized regions considered in section II. Nearly all 
the radiation beyond the Lyman limit, and also that in the Lyman 
lines L6, Lv, etc., has here been converted into La and low-frequency 
lines. Free electrons, supplied by elements of lower ionization po- 
tential than hydrogen, will all have relatively low kinetic energy. 
Therefore, in these parts of interstellar space there is almost no 
excitation of the state m = 3 or of higher states. The state n = 2 
is probably superexcited by a factor of 10, very roughly, owing to the 
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conversion of continuous radiation into La, while the metastability 
of the substate 2s is of no consequence here. 

Eddington’ has suggested that hydrogen in un-ionized regions of 
interstellar space is present in molecular form. It is not intended to 
attempt a complete solution of this problem here, but a few points 
may be noted. If the equilibrium between hydrogen atoms and mole- 
cules is determined principally by one process and its reverse proc- 
ess, the ratio of atoms to molecules will be given by an equation 
analogous to (1). (In interstellar space the lowest vibrational and 
rotational level only will be excited. The reduced mass 3my occurs 
instead of m., however.) If the frequency of the radiation absorbed 
in the dissociation process considered were less than that of the Ly- 
man limit, the factor corresponding to e~™« would be of the order of 
magnitude 1. In that case the ratio of molecules to atoms would be, 
roughly, 10-°.V; i.e., the number of molecules would be negligibly 
small. If, however, the frequency of the dissociation process is great- 
er than that of the Lyman limit—and this is certainly the case if the 
dominant dissociation and recombination processes are the same as 
in the laboratory—then the number of dissociations is practically 
zero. On the other hand, the reverse recombination process, which 
involves the collision of one normal and one excited hydrogen atom, 
is extremely slow in interstellar space. Using the approximate value 
for the effective cross-section of this process given by Terenin and 
Prileshajewa™’ and by Rompe," namely, 10-5 times the gas-kinetic 
cross-section, it is found that every hydrogen atom collides to form 
a molecule once in (very roughly) 10*t years. This means that mole- 
cule formation during any accepted time-scale is negligible. This, 
of course, does not definitely settle the question. It should be added 
that the mixture of material in the galactic system will tend ulti- 
mately to bring any molecules that may exist into ionized regions, 
where they will quickly dissociate. A sufficiently strong mixture will 
thus lead to practically complete dissociation of molecules into 
atoms. 

IV 

According to the analysis of the preceding sections, one would ex- 

pect—if the density of interstellar hydrogen is within a certain range 


13 Phys. Zs. Soviet Union, 3, 337, 1932. ™4 Phys. Zs., 37, 807, 1936. 
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—to find Balmer-line emission limited to certain rather sharply 
bounded regions in space surrounding high-temperature stars or 
clusters of high-temperature stars, the emission being of uniform 
strength, without any concentration toward particular stars, within 
these regions. With a view to the relative ionizing efficiency of O 
and B stars, as discussed in section II, one would identify the high- 
temperature stars in question principally with O stars. 

Comparing this picture with the observations by Struve and 
Elvey,' referred to at the beginning of this article, of extended re- 
gions of Ha emission, it is tempting to interpret the observed areas 
of uniform Ha emission as the projections of bounded regions of in- 
terstellar space ionized principally by clusters, or condensations, of 
O stars. We shall, for the present, accept this interpretation as a 
working hypothesis. 

Our working hypothesis may be regarded as a specialized version 
of the more general hypothesis put forward by Struve? in the investi- 
gation previously referred to. 

Several clusters or clusterings of O stars are known.” In accord- 
ance with our working hypothesis we connect the observed extended 
Ha emission region in Cygnus with the cluster of O stars having its 
center at /, = 42°, b, = +1°. About ten of the known O stars pre- 
sumably belong to this cluster. Likewise, we connect the extended 
Ha emission in Cepheus with the condensation of O stars around 
galactic longitude 70°. 

We can now immediately make use of the analysis of section II] 
to derive the density of interstellar hydrogen in the regions in ques- 
tion. Struve’?, from the observed strength of the Ha emission in 
Cygnus and Cepheus, derived V, = 3 - 10°" cm~', assuming the 
effective length of the emitting column to be tooo parsecs. Using 
this number we estimate, with the aid of Table 6, that the density 
N of hydrogen is equal to 2-3 cm~3. We have here assumed that 
practically all the free electrons result from hydrogen (V”" = N.,, 
and also N”’ = N, the ionization being nearly complete). This 
is probable a priori and is further supported by reasons stated be- 
low. With an estimated depth of the emitting column of 300 parsecs, 

15 C, H. Payne, ‘‘The Stars of High Luminosity,”’ Harvard Obs. Monograph, No. 3, 
chap. vi, Cambridge, Mass., 1930; J. S. Plaskett, Pub. Dom. Ap. Obs., Victoria, 2, 339, 
1924; O. Struve, A.N., 231%, 17, 1027. 
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instead of 1000 parsecs, one would find NV equal to about 4 cm-3. The 
resulting V is proportional to the square root of the assumed value 
of NV. 

Assuming, for the present, that NV = 3, we find, with the aid of 
Table 5, that the diameter of the region ionized by a cluster of n 
stars of the average spectral type O7 is about 80 n'’3 parsecs. With 
n = 10 (cf. above) we find a diameter of about 200 parsecs. This 
agrees about as well as could be expected with the estimated diam- 
eter of the Cygnus emission regions of, say, 200-300 parsecs. 

The total intensity of a hydrogen-line emission area is propor- 
tional to the number of ionizing O stars. The total area is propor- 
tional to n?’3, the depth, and hence the surface intensity, to n'/3. 

With regard to the extension of the areas of Ha emission in the 
direction at right angles to the galactic plane, consideration of the 
concentration of hydrogen toward the galactic plane is of impor- 
tance. Calculations based on the assumption of a static interstellar 
hydrogen cloud may perhaps serve as a rough guide. Using the em- 
pirical field of force in the direction at right angles to the galactic 
plane as derived by Oort," it is found that the density of hydrogen 
decreases with the height above the galactic plane approximately ac- 
cording to a Gaussian error law, with a dispersion of about 1.3 
parsecs X 7*/?,where 7 is the temperature of the hydrogen gas. With 
T = 10,000°, a dispersion of 130 parsecs results. Though the true 
distribution will probably differ from this, owing to the effects of 
turbulence, the order of magnitude should be correct. 

It is, of course, very important to find how large a fraction of the 
relevant galactic stratum of hydrogen is ionized. From the esti- 
mated density of O and B stars, in connection with Table 5, assum- 
ing V = 3 cm-, we are led to the very rough estimate that about 
one-tenth of the stratum is ionized. The numerical value of the 
fraction is inversely proportional to the square of the density of in- 
terstellar hydrogen. 

As the observational material on extended areas of hydrogen-line 
emission accumulates, it will be possible to test the working hypothe- 
sis formulated in this section. The necessity of a critical test is ob- 
vious, particularly in view of the fact that a revision of the tempera- 
ture scale might lead to a somewhat different picture of the relative 


%© Bull. Astr. Inst. Netherlands, 6, 249, 1932. 
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importance of O and B stars and that a decrease in the accepted value 
of N might lead to a less pronounced separation between ionized and 
nonionized regions. 

Further observations of hydrogen-line emission areas would be 
especially valuable, as they would lead to an improved knowledge 
of the relative extent of that part of the galactic stratum in which 
hydrogen is ionized. 

At the present moment it can be said that a number of determina- 
tions, or estimates, of the density of interstellar hydrogen lead to 
values in approximate mutual agreement. First, Dunham’s'’ dis- 
covery of interstellar lines of neutral calcium and the discussion of 
this phenomenon by Dunham and Struve suggest as the most plausi- 
ble choice of the value of the hydrogen density the maximum value 
compatible with the results of Oort’s dynamical investigations, i.e., 
N = 2-4 cm™3 (cf. Struve?). The analysis based upon the strength 
of interstellar Ha emission leads to a value of NV around 3 cm~. The 
observed extent of Ha emission areas confirms the approximate cor- 
rectness of this value. Finally, it may be mentioned that Edding- 
ton’s estimate, based upon the average separation of nebular con- 
densations, when modified to correspond to a mean molecular weight 
of 1, leads to N = 5 cm". 

V 

We shall consider briefly the state of ionization of elements other 
than hydrogen. In those‘regions of interstellar space where hydrogen 
is ionized, the electron density V, is equal to N. All ionizing radia- 
tion of wave lengths shorter than that of the Lyman limit comes 
from stars within the ionized region, i.e., particularly from the same 
stars which ionize the hydrogen. This radiation is not very much 
reduced by hydrogen absorption except near the boundary of the 
ionized region (cf. sec. II). In the case of sufficiently abundant ele- 
ments of high ionization potential the ionization may be checked in 
a way analogous to that in which the ionization of hydrogen is 
checked. With obvious modification, the analysis in section II holds 
for such cases. Thus, it is seen that, if helium were as abundant as, 
or more abundant than, hydrogen, ionized helium would be restricted 
to a volume considerably smaller than that of ionized hydrogen. If, 
on the other hand, helium is considerably less abundant than hydro- 


17 Pub. A.S.P., 49, 26, 1937; Nature, 139, 246, 1937. 
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gen, the ionization of He1to Hew will usually be almost com- 
plete up to the boundary of the region of hydrogen ionization; but in 
that case, of course, helium contributes considerably less to the num- 
ber of free electrons than hydrogen. This argument tends to con- 
firm the validity of our assumption that the free electrons are con- 
tributed mainly by hydrogen in those regions where hydrogen is 
ionized. It also shows that helium D,-line emission is handicapped 
in comparison with hydrogen Balmer-line emission. 

Struve and Elvey' observed the emission line \ 3727 of [O 1] in 
many of the areas where Ha was seen, while NV, and N, of [O 111] were 
rarely observed. As the mechanism of excitation of the lines is 
similar, this probably indicates that oxygen is present mainly as 
Ou. Further observations on the occurrence of V, and N, may lead 
to valuable results, particularly so if it should turn out that the abun- 
dance of oxygen is high enough to restrict the ionization of O 11 to 
O 111 (ionization potential 34.9 volts) by a mechanism similar to the 
hydrogen mechanism. 

Ultimately, the observed relative intensity of interstellar hydro- 
gen and oxygen emission lines will lead to a reliable value of the rela- 
tive abundance of interstellar hydrogen and oxygen. A very rough 
estimate is easily made. The effective cross-section for the excita- 
tion of the upper state of \ 3727 by electron collision is estimated to 
be somewhat smaller than 10~'? cm’. Assuming that somewhat less 
than one-tenth of the electrons have sufficient energy to excite 
d 3727, and putting the cross-section for capture of a free electron 
to the state m = 3 equal to about 10~*' cm’, we derive from the ob- 
served fact that \ 3727 is somewhat fainter than Ha that there are, 
roughly, 107-103 H ions to every Otrion. Since hydrogen is largely 
ionized, and oxygen is mostly present as O 1, the resulting inter- 
stellar abundance is 10-?-10 * oxygen atoms to every hydrogen atom. 

We consider next the ionization of elements with ionization po- 
tential less than that of hydrogen, in those regions where hydrogen 
is ionized. In this case the ionizing radiation comes from the entire 
galactic system, so that the customary calculations of the dilution 
factor and the effective temperature hold, while the electron density 
Nis equal to NV. It should be noted that the O stars which ionize the 
hydrogen do not contribute especially to the ionization of such ele- 
ments as calcium and sodium. 
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Outside the regions in which hydrogen is ionized, there will be 
almost no ionization for ionization potentials larger than that of 
hydrogen. Of the most abundant elements, helium, nitrogen, and 
oxygen will be un-ionized, while carbon and the metals will be ion- 
ized. Of the abundant metals, only calcium will be doubly ionized. 
It may be noticed that titanium will be present mostly as 77 1. This 
is of interest in connection with Dunham’s” observations of interstel- 
lar 77 11 lines. 

It is quite likely that in the regions of un-ionized hydrogen most 
of the free electrons are furnished by carbon. With an interstellar 
carbon abundance of the same order as the abundance derived for 
oxygen, we find that NV, is about 107103 times smaller than in the 
regions where hydrogen is ionized. This figure is, however, as yet 
extremely uncertain. It is probably best to postpone the detailed 
discussion of the ionization of the metals until two basic constants 
are known with greater certainty, namely, the fraction of the galactic 
stratum in which hydrogen is ionized, and the ratio of electron den- 
sity in regions of ionized and un-ionized hydrogen, respectively. With 
the present values for these constants it is found that interstellar 
Ca 1 lines are almost entirely absorbed in the regions of ionized hydro- 
gen, while Ca 11, Na 1, and K 1 lines are also mostly absorbed here. 

Restriction of the absorption of an interstellar line to a fraction 
of the galactic stratum naturally affects the abundances calculated 
from the observed intensities of the line. Increasing the abundance 
of calcium and sodium found by Struve? (forming an average abun- 
dance of the two elements) on this account, and using V, = 3 in the 
absorption region, we arrive at a density of 3.10->cm~. Adopting, for 
interstellar space, the relative cosmic abundances of oxygen, calcium, 
and sodium, as determined by Goldschmidt,"* the corresponding 
density of oxygen is found to be 5.10-? cm~3. The resulting inter- 
stellar abundance is 2 - 10~? oxygen atoms to every hydrogen atom, 
which is of the same order as the value determined above from the 
strength of the line \ 3727. It is hardly necessary to emphasize the 
fact that some of the factors entering into this calculation will 
certainly need considerable revision.’ 

18 Norske Videnskaps-Akademi, Oslo, Mat.-Naturv. KI., No. 4, 1937. 


19 In a more refined calculation of the ionization of an element in interstellar space 
the effect of recombination processes to excited states, followed practically always by 
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Struve” has recently emphasized the relatively much greater 
sensitivity of the interstellar absorption-line method as compared 
with the interstellar emission-line method. Even so, it seems already 
quite clear that the emission-line method is capable of yielding, 
ultimately, extremely valuable results with regard to the determina- 
tion of the interstellar density of the most abundant elements, all of 
which show no interstellar absorption lines in the accessible wave- 
length region. 


The author is much indebted to Dr. O. Struve for placing at his dis- 
posal the results of investigations on the subject of this paper before 
their publication. 


COPENHAGEN OBSERVATORY 
February 15, 1939 


cascading to the ground state, must be taken into account. The effect will always tend 
to decrease the degree of ionization. For hydrogen the effect is not very marked, leading 
to an extra factor on the right-hand side of equation (1) estimated to be about equal 
to 3. This would correspond to a decrease of the so-values calculated in the present pa- 
per by about 30 per cent. For elements like sodium and potassium, however, the 
effect is very marked. It is well known (cf. e.g., n. 4) that for these elements the ab- 
sorption from the ground state is much smaller than for hydrogen, while, according 
to an unpublished investigation by Mr. M. Rudkjébing, the excited states are approxi- 
mately hydrogen-like with regard to continuous absorption. It follows that for these 
elements the number of recombinations to excited states is far greater than the number 
of recombinations to the ground state. In interstellar space the number of dissociation 
processes is, of course, proportional to the continuous absorption coefficient for the 
ground state. We conclude that the degree of ionization of sodium and potassium in 
interstellar space is quite considerably reduced by the effect in question. The abun- 
dances of sodium and potassium calculated from the strength of interstellar absorption 
lines are similarly reduced. According to an unpublished result by Dr. L. C. Green, 
kindly communicated to the writer, Ca 11 is intermediate between sodium and hydrogen 
with regard to the continuous absorption from the ground state. Hence, a reduction of 
the calculated calcium abundance is indicated, but it is smaller than for sodium. Pend- 
ing more accurate calculations, it is suggested that the remainder (cf. Struve?) of the 
well-known discrepancy with regard to the relative strength of interstellar Va 1 and Ca 11 
lines will disappear when the effect just considered is taken into account. For Ca 1 the 
continuous absorption coefficient from the ground state has not yet been calculated. 
It is very probable, however, that it is considerably smaller than for hydrogen, the 
excited states being more hydrogen-like with respect to continuous absorption. If 
this is the case, it follows that the interstellar ratio of Ca 1 to Ca 11 is considerably larger 
than that calculated from the customary equation. This would tend to reconcile the 
value of the electron pressure derived from the observed strength of Ca1 and Ca 11 (cf. 
Dunham"? and Struve?) with that found in the present investigation. 


Zs. f. Ap., 17, 316, 1939. 
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TWO NEW WHITE DWARFS; NOTES ON 
PROPER MOTION STARS 
ABSTRACT 

This note is a brief report on some of the observations made with the 82-inch tele- 
scope during the first month of its operation. 

During the first month of regular operation the 82-inch telescope 
was used mainly in the determination of accurate spectral types for 
stars of large proper motion (not less than 0% 30 per year). About 400 
spectra were obtained for about 250 stars between visual magnitudes 
8 and 16. Three classes of objects of special interest are expected to 
be found in such a survey: (1) white dwarfs; (2) intermediate white 
dwarfs or, more generally, stars not over 2 or 3 mag. below the main 
sequence; and (3) stars of large (spectral) parallax. It is found that 
the second group extends almost along the whole main sequence. 
Since these stars merge into the main sequence and are much more 
similar to main-sequence stars than to white dwarfs (probably also 
in the interior),’ the name ‘“‘subdwarfs”’ is suggested for this class 
of stars, in analogy with “subgiants.”” This name will prevent the 
confusion of these stars with the white dwarfs proper which are very 
much fainter. 

During the past month two new white dwarfs were found, perhaps 
two dozen of subdwarfs, and more than a dozen stars with spectral 
parallax in excess of o71. These incidental results are not regarded 
as the sole object of the survey; its main purpose is the derivation 
of the properties of the near-by stars by completely covering the 
known stars of proper motion above 0730. About 600 faint stars 
have been observed so far since the work was started last year at the 
Yerkes Observatory. 

Table 1 contains the white dwarfs and the most probable sub- 
dwarfs found with the 82-inch telescope. In addition, most of the 
known white dwarfs were observed (a few objects were in inaccessi- 
ble parts of the sky) and a summary of the comparisons between 
these objects is given below. 

Wolf 457 has a color index of about zero, and has a continuous 

t This may be seen, e.g., from Table 11, p. 508, of B. Strémgren, Erg. d. Ex. Nat., 16, 


1937: 
548 

















NOTES 549 


spectrum between about 6600 and 3800 A. Ross 808 has strong and 
broad Balmer lines and apparently belongs to that group of the 
white dwarfs of which 4o Eridani B is the prototype. According to 
the writer’s observations, 8 of the white dwarfs now known belong 
to that group. Their spectra are almost indistinguishable and show 


TABLE 1 





| | 
Star a 5 Myis P.M. | Sp 1oo7(Sp) 
Wolf 457 $2266 23 + 4:02 15.6 1705 O 
Ross 808 ESS). 7 +37 06 14.4 0.58 \o 
Ross 338 2 51.2 +48 28 13.2 0.42 G2 0” 20 
Ross 585 o 94A +33 43 12.9 1.41 K4 0.83 
Ross 600 4 35.4 +22 43 13.0 0.68 Ko 0.38 
Sel. A. 342 5 60.7 +60 2 14.8 0.50 K3 0. 26 
Ross 40 5. 30-4 + 9 13 11.9 0.60 F8 0.29 
Ross 988 ?. 25. (5 +38 II 13.8 0.74 Ko 0.28 
Ross 390 7 20.1 —10 10 11.9 0.62 Gs 0.42 
Wolf 1059 at +29 18 12.6 0.53 Go 0.2 
AC 78°3159 8 40 +78 54 12.3 0.61 G8 0.46 
Ross 890 9 44.5 + 7 05 12.4 ©. 32 F8 C2) 
Ross 892 [oO 22-3 + 1 55 II .4 0.38 F8 0°. 36 
Ross 108 [1 00.3 +53 45 Es. 2 O: 77 K6 0.46 
Sel. A. 33 II 30.9 +29 24 13.8 0.88 <3 0.42 
Ross 114 [0-308 +14 57 13.8 a G5 0.17 
Ross go4 iI 36.1 +25 39 12.9 0.47 G8+ 0.35 
+26°2251 II 30.5 +26 07 10:6 | oO. 52 F8 0.52 
+51°1696 II 41.4 +51 27 g.§ | 1.02 Go 1.0 
Ross go8 II 47.5 +28 06 2.1 0.92 G5 0.38 
—4°3208 12 02.2 — 5 10 10.2 ©.54 Fo 0.35 
Rose 453....| %2 05.8 +o 58 EE.2 0.50 Fot °0.22+ 
2°2538 12 20.5 + 15! 9.8 0.51 F2 0.48 
Ross 974 12 56.9 — I 33 12.5 0.57 G5 ©.32 
Ross 838 13 56.8 + 9 25 ee 0.88 G5 0.50 
Ross 841 rs $7.3 — 5 10 11.2 0.44 A2s 0.13 























only the strong Balmer lines up to about H¢. These stars are, in 
order of right ascension: 1166 h Persei,? 40 Eridani B,’ Hertzsprung 
3=CC 308,! Furuhjelm I, 182,5 Wolf 485 = BD —7°3632,° Ross 808,’ 

2 Oosterhoff, B.A.N., 6, 39, 1930. 3 Adams, Pub. A.S.P., 26, 198, 1914. 

4 Kuiper, Pub. A.S.P., 49, 341, 1937; see also Ramberg, A.N., 265, 111, 1938. 

5 Deutsch, Observatory, 62, 53, 1939. The proper motion was known; the spectral 
class used was from the Bergedorf Spektral Durchmusterung. 


6 Kuiper, op. cit., 47, 307, 1935. Dr. A. van Maanen has pointed out to the writer 
that the co-ordinates of BD —7°3632 agree with those of Wolf 485, so that the BD desig- 
nation may be used in spite of the faintness of the star, 12.0 mag. 


7 This note. 
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Wolf 1346,° and AC 82°3818.9 In addition, the two white dwarfs in 
the Hyades'® probably belong to this group. The absolute magni- 
tude of these objects seems to be about +11; a recent parallax de- 
termination at the Cape of Wolf 485 makes the weighted mean paral- 
lax of this star 0%064+8, giving an absolute magnitude of +11.0, 
nearly the same as for 40 Eridani B, +11.2. 

Ross 627" does show the Balmer lines, but they are much shallow- 
er than for the stars of the preceding group. This fact may be in- 
terpreted as being due to a larger surface gravity—a result of the 
lower absolute magnitude (133-14) and probably a much smaller 
radius. This star forms a transition to the stars with practically 
continuous spectra, in which the Balmer lines have become so shal- 
low as to be invisible. The three cases in this group, all found by the 
writer, are AC +70°8247," Wolf 219,8 and Wolf 457.7 Shallow lines 
may exist in these spectra, but it is obvious that these spectra are 
unlike any found for main-sequence stars. The fact that the Balmer 
lines are invisible can be interpreted as the result of high surface 
gravity,” possibly also by rapid rotation. The ad hoc assumption 
that no hydrogen exists in these stars, even in the atmosphere, which 
as a result of excessive surface gravity may contain a much greater 
abundance than the star as a whole, seems very unsatisfactory in 
view of the great strength of the Balmer lines in the majority of the 
white dwarfs. (The theoretical difficulty with the energy genera- 
tion is not serious unless ‘it is shown that there is a thorough mixing 
of the elements throughout the star; on the contrary, this very diffi- 
culty points to the opposite in connection with the strong Balmer 
lines in the majority of the white dwarfs.) 

In addition, five other white dwarfs are known: the companion to 
Sirius (A5),"4 that to Procyon (spectrum unknown),'’ van Maanen 

8 Kuiper, of. cit., 46, 290, 1934. 9 Tbid., 47, 279, 1935. 

% One was suspected by Schwassmann and van Rhijn, M.N., 94, 517, 1934. This was 
confirmed and another added by Ramberg in A.N., 265, 111, 1938. 

« Humason and van Maanen, Pub. A.S.P., 48, 179, 1936. 


2 Op. cit., 47, 307, 1935- 

13 [bid., p. 96, 1935. The writer has not yet been able to observe this star with the 
82-inch telescope, but Mr. Humason, of the Mount Wilson Observatory, has informed 
him that he found Wolf 219 to have a continuous spectrum. 


14 Adams, Pub. A.S.P., 27, 236, 1915. 5 E.g., Ap. J., 88, 490, 1938. 



































NOTES 551 


2 (Go; actually the spectrum differs from any type applicable to 
ordinary dwarfs),"®° Wolf 489,'7 and Ross 451."7 

Wolf 489, with absolute magnitude about 15.1 visually, 15.8 pho- 
tographically, as derived at Mount Wilson, is perhaps the most pe- 
culiar of the white dwarfs. At the Yerkes Observatory its spectrum 
was found to be earlier than M; later (unpublished) it was found 
not to be of type K. Two spectra obtained with the 82-inch tele- 
scope, covering the region between 3800 and 6600 A, fail to show 
any lines, although the intensity distribution is clearly in agreement 
with the approximate value of the color index derived (about 0.7 
mag.). The most amazing fact is the absence of the H and K lines 
which would be expected to be the strongest lines in the spectrum. 
The presence of faint traces of these (and other) lines is, of course, 
not precluded, but it is certain that none of the strong features used 
in the Draper classification or used for classification in the visual 
region is present in anything approaching their normal strength. 
The color index of the star will be determined more precisely in order 
to obtain a good determination of the radius of the star. 

Ross 451 (K4) is thus the latest among the white dwarfs; but 
since it seems to be only between 3 and 4 mag. below the main se- 
quence, it might perhaps more properly be considered a subdwarf. 
Further determinations of its parallax will be important in this con- 
nection. 

The probable subdwarfs in Table 1 are selected mainly on the 
basis of abnormally high transverse motion computed under the as- 
sumption of the star being a normal dwarf. If the proper motion is 
one hundred times the spectral parallax or more, the star is probably 
below the main sequence. In some cases the ratio is about 300; then, 
even if the transverse motion is as high as 474 km/sec, the star must 
be over 2 mag. below the main sequence. Actually small observa- 
tional errors will make the distinctions less precise, and some of the 
stars can only be presented as probable cases. There is, however, an 
additional criterion for identifying these subdwarfs for types earlier 
than about G5: the weakness of the Balmer lines and the G band. 
A more extensive discussion of these stars will be published later. 

One new double star was found among the near-by stars: +10° 


© Seares, Pub. A.S.P., 30, 191, 1918. 17 Kuiper, Ap. J., 87, 593, 1938. 
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1857, 975, K6+; separation 2”, Am = 3.5. The star has a third 


component, separation 115’’, 122, M4+. The spectral parallax of 
the system is 0708. 
G. P. KUIPER 
McDOoNnaLpD OBSERVATORY 
April 1939 





THE BEHAVIOR OF He 15015 IN 
¢ TAURI AND @¢ PERSEI 
ABSTRACT 

The Het line \ sors (2'S — 3'P), predicted to be sharp in ¢ Tauri and ¢ Persei, as 
a result of dilute radiation in the nonrotating outer shells of these stars and the metasta- 
bility of the 2'S level, is found to be as broad as D3 in ¢ Tauri and absent or present 
only as a very faint, broad line in @ Persei. 

In a recent paper on excitation phenomena in outer shells of 
stars Struve and Wurm’ explain the observed strength and sharpness 
of the Het line \ 3965 (2'S — 4'P) in ¢ Tauri (B3e) and ¢ Persei 
(Bone) as being due to the metastability of the 2'S level and to 
the presence of nonrotating outer shells. Dilute radiation presum- 
ably causes the 2'S level to become richly populated with respect 
to the higher nonmetastable levels, thus accounting for the strength 
of the line, while its sharpness is considered a result of the nonrotat- 
ing shell. The corresponding lines arising from the photosphere, pre- 
sumed to be rotating rapidly, are greatly broadened and indistinct 
but bear the expected relative intensities as a consequence of the 
absence of dilute radiation. 

The only other available line arising from the 2'S level is \ 5015 
(2'S — 3'P). It was not observed by Struve and Wurm, owing to its 
presence in the green part of the spectrum, but it was assumed that 
it would behave similarly to » 3965. 

Two plates in the visual region of ¢ Tauri (October 5 and Decem- 
ber 24, 1936) and one of @ Persei (September 22, 1937), taken by 
Dr. Cherrington with the Yerkes auto-collimating grating spectro- 
graph (27 A/mm), were available in the files of spectrograms at the 
Perkins Observatory. Plate XX V shows X 5015 and D3 in (a) ¢ Tau- 
ri and (b) @ Persei. Microphotometer tracings of these regions are 
shown immediately below. 

Inspection shows that in ¢ Tauri \ 5015 appears as broad as D3 


1 Ap. J., 88, 84, 1938. 
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PLATE XXVI 


THE BriIGHT NEBULOSITY SOUTH OF S MONOCEROTIS 
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of Het and that \ 4922 and X 5047 are present and have similar 
profiles to that of A 5015. The singlet He 1 \ 6678, not shown, hav- 
ing the same laboratory intensity as \ 5015, appears in the star with 
much the same strength and profile as \ 5015. It is possible that a 
narrow faint absorption line is superposed over the broad line. If 
this line is real it is probably due to the shell, but the evidence is not 
complete. 

Struve and Wurm report that in @ Persei 3965 is always sharp. 
We should expect, therefore, that \ 5015 should also be sharp in this 
star. Actually, however, it is absent or present only as a very faint, 
broad line. The other singlet lines in the visual region (AA 4922, 
5047, and 6678) likewise appear absent or, at most, very faint. The 
triplet D3 line, however, is fairly strong. This weakness of the sin- 
glets relative to the triplets in ¢ Persei is confirmed, but the fact that 
\ 5015 behaves as do the other singlet lines is significant and puz- 
zling. The emission line near \ sors in ¢ Persei is \ 5018 of ionized 
iron, and it may affect the appearance of \ 5015. The lines A 4923 
and \ 5169 (a°S — z°P) of Fe 11 are also present in emission. It may 
be noted in passing that these lines have about the same relative in- 
tensities in @ Persei as they had in Nova Herculis,? \ 5169 being dis- 
tinctly the strongest. 

In absorption these lines are reported about equal by Marshall? 
in the late B- and early A-type stars but weak or absent in the early 
B stars. In ¢ Tauri \ 5169 is very weak, but there is a distinct emis- 
sion at \ 5018 which may affect the appearance of \ 5015 by blend- 
ing. However, the evidence suggests that \ 5015 is not appreciably 
different in behavior from the other singlet lines of helium in ¢ Tauri 
and @ Persei. The behavior of \ 3965, arising from the same level, 
becomes thus still more puzzling. J. A. Hynex 


PERKINS OBSERVATORY 
April 1939 





THE NEBULOSITY NEAR S MONOCEROTIS 


The photograph reproduced in Plate XXVI was obtained on 
February 15, 1939, with the 24-inch reflector of the Yerkes Observa- 
tory. The emulsion used was Eastman Ortho Press. The central 


2 Contr. Perkins Obs., No. 2, p. 24, 1935. 3 Ap. J., 82, 97, 1935. 
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portion of the nebula shown in Dr. Van Biesbroeck’s photograph, 
Plate XXVII, is overexposed. An interesting feature of this nebula is 
the comet-shaped dark marking south of the main nebulosity which 
shows a conspicuous luminous rim of the type described by Struve.' 
W. W. Morcan 
FRANCES SHERMAN 
YERKES OBSERVATORY 





THE NEBULOSITY NEAR S MONOCEROTIS 
In the course of the first tests of the 82-inch telescope the nebulo- 
sity near S Monocerotis (P]. XX VII) was photographed at the prime 
focus on March 13, 1939, with an exposure of five minutes. This bright 
nebula is particularly interesting because of the conspicuous absence 
of nebulous matter near two of the stars imbedded in the nebula. 
This is a well-known phenomenon in several other emission nebulae. 
G. VAN BIESBROECK 
McDONALD OBSERVATORY 


PHOTOGRAPH OF A DIFFUSE NEBULA IN ORION 

A small nebula, photographed by H. A. Lower' has been examined 
by Struve and Elvey? with the nebular spectrograph of the McDon- 
ald Observatory (No. 43 of their list). It has been photographed at 
the Yerkes Observatory with the 24-inch reflector on Agfa Super- 
Pan Press film. A Wratten Ciné-Red filter was used to cut out 
radiation from the night sky to the violet of \ 6000. An exposure of 
two hours on March 17, 1939, shows the lumpy central part, which, 
as Struve and Elvey point out, is considerably brighter than the rest 
(Pl. XXVIII). A faint extension northward is also shown. The star 
BD-+15°1079 lies within the brighter portion; BD+15°1075 lies a 


little outside. 
Joun A. O'KEEFE 
YERKES OBSERVATORY 
March 1939 





ERRATUM 
In the January issue, page 121, in Table 1, No. 28, under “Identification”’: 
Ross 26 should read Ross 25. 


1 Ap. J., 85, 208, 1937. 
1 Ap. J., 89, Pl. XIV, 1939. 2 Ap. J., 89, 1939. 














PLATE XXVII 


THE BRIGHT NEBULOSITY SOUTH OF S MONOCEROTIS 


























PLATE XXVIII 











Lower’s DIFFUSE NEBULA IN ORION. SCALE 20” PER MILLIMETER 

















